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Context. Observations of the formation and evolution of massive galaxy clusters and their matter components provide crucial con- 
straints on cosmic structure formation, the thermal history of the intracluster medium (ICM), galaxy evolution, transformation pro- 
q , cesses, and gravitational and hydrodynamic interaction physics of the subcomponents. 

5_l ■ Aims. We characterize the global multi- wavelength properties of the X-ray selected galaxy cluster XMMU J 1230.3+ 1339 at z= 0.975, 

^ ' a new system discovered within the XMM-Newton Distant Cluster Project (XDCP). We measure and compare various widely used 

| mass proxies and identify multiple cluster-associated components from the inner core region out to the large-scale structure environ- 

ment. 

Methods. We present a comprehensive galaxy cluster study based on a joint analysis of X-ray data, optical imaging and spectroscopy 
observations, weak lensing results, and radio properties for achieving a detailed multi-component view on a system at z ~ 1. 
k>" ' Results. We find an optically very rich and massive system with M 20 o - (4.2 ± 0.8) x 10 14 Mq, r x2 5oo - 5.3^g keV, and 

_Lx°500 (6.5 ± 0.7) x 10 44 ergs~'. We have identified a central fly-through group close to core passage and find marginally extended 
1.4 GHz radio emission possibly associated with the turbulent wake region of the merging event. On the cluster outskirts we see 
evidence for an on-axis infalling group with a second Brightest Cluster Galaxy (BCG) and indications for an additional off-axis group 
accretion event. We trace two galaxy filaments beyond the nominal cluster radius and provide a tentative reconstruction of the 3D- 
. accretion geometry of the system. 

1 Conclusions. In terms of total mass, ICM structure, optical richness, and the presence of two dominant BCG-type galaxies, the newly 

' confirmed cluster XMMU J1230. 3+1339 is likely the progenitor of a system very similar to the local Coma cluster, differing by 

7.6 Gyr of structure evolution. This new system is an ideally suited astrophysical model laboratory for in-depth follow-up studies on 
the aggregation of baryons in the cold and hot phases. 



X 



Key words, galaxies: clusters: general - galaxies: clusters: individual: XMMU J1230.3+1339 - X-rays: galaxies: clusters - galax- 
ies:ellipticals and lenticular - galaxies: evolution - cosmology: dark matter - cosmology: observations 

1 . Introduction evolution of Dark Energy (DE). This potential has recently been 

demonstrated based o n moderate size X-ray cluster samples ou t 

Massive galaxy clusters at high redshift are unique laboratories tQ redshift uni ( [yi khlinin et al . 2 009fi lAllen et al J 12003) 

to study galaxy evolution in the densest environments, the ther- and k reflected m the ^ m ^ a of ongoing or planned galaxy clus . 

modynamic properties and chemical enrichment of the hot intra ter surveys in var i ous wavelength regimes, 
cluster medium (ICM) at large look-back times, and the structure 

and growth of the underlying virialized Dark Matter (DM) halos. Progresses in observational capabilities and search tech- 

Moreover, distant galaxy clusters are among the most promis- mc l ues over the P ast five y ears have led to an increase 

ing cosmological probes to shed new light on the properties and ln the number of spectroscopically confirmed z £ 1 clus- 



* Based on observations obtained with ESO Telescopes at the Paranal tems - ° n me X " ra y side > new identifications of X-ray lumi 
Observatory under program ID 078.A-0265 and 081.A-0312, and ob- nous clusters are currently largely driven by XMM-Newton 



servations with the Large Binocular Telescope (LBT), and the X-ray serendipitous surveys (e.g. iMullis et al 
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1201 Oh an d smaller dedicated surveys of a few square de- 
grees ( e.g. lAndreon et af]|2005b [Bremer et aDl2006t iPierre et al.l 
120061: iFinoguenov et al.l l2007i 120101: iTanaka et al.l |2o7o1k 
Optical search methods b ased on overdensities of galax- 
ies at similar color (e.g. G ladders & Yed 120051) have re- 
cently been successfully extended to t he near-infrared (NIR) 

and the mid-infrared 
Eisenhardt et all 



regime (e.g. Ivan Breukelen et al. 1 120061) 
(MIR ) with Spitzer (e.g. [Stanford et all 
20081: IWilson et al.l 120091: iMuzzin et al 
20101) . The first detections of ga laxy 
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Demarco et al.l 
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the Sunyaev-Zeldovich effect (SZE) (IStaniszewski et all 120091: 
IVanderlinde et al.l2010l) have opened another promising window 
for the selection and study of massive distant systems. 

A crucial pre-requisite for linking galaxy cluster survey re- 
sults to cosmological model predictions is reliable and well 
quantified mass-observable relations. In this respect, weak grav- 
itational lensing techniques have emerged as a promising tool to 
achieve a robust calibration of various scaling relations. In the 
low redshift Universe, weak lensing (WL) based mass measure- 
ments applied to sizable cluster samples have established mass- 
observable relations with significantly improved accurac y (e.g. 
iReves et al.ll2008T: IZhang etaT]|2008b iMarrone et al.ll2009b . First 
weak lensing constraints on the evol ution of scaling relation s for 
X-ray groups are now available (e.g. lLeaufhaud et al.l2010t) . and 
significant WL signals have b een reported for indiv idual mas- 
sive clusters out to z ~ 1.4 (e.g. lJee et alj|2006ll2009l) . However, 
all mass estimation metho ds are potentially influ enced by dif- 
ferent inherent biases (e.g. lMeneghetti et al.ll2.Q10l) emphasizing 
the importance of an inter comparison of different techniques, in 
particular for high redshift systems. 

Significant progress is yet to be made to establish robust 
mass-observable relations at z ~ 1 and beyond and to compare 
vices and virtues of mass proxies derived from optical/IR, X-ray, 
WL, and SZE observations at these redshifts. The current limita- 
tion is the still persisting small number of known galaxy clusters 
at z > 0.9 that are sufficiently massive to enable the detection and 
cross-comparison of all observational techniques on an individ- 
ual cluster basis. Using concordance model parameters, we can 
estimate the surface density of suitable test objects at z > 0.9 
with M200 ^4xl0 14 M(7) to be abou t one object in a solid angle 
of 100 deg 2 (see e.g. lJee et al.ll2009l) . 

In this paper we present first details of the newly identified 
cluster of galaxies XMMU J1230.3+1339 at redshift z = 0.975. 
We perform a cluster characterization based on a joint analysis 
in X-rays, optical imaging and spectroscopy, and weak lensing. 
The three main science objectives of this work are (i) a multi- 
wavelength characterization of the physical system parameters, 
(ii) the derivation and cross-comparison of a dozen different 
total mass proxies, and (iii) the identification of substructure- 
and large-scale structure components for an assessment of the 
dynamical stat e of th e system. An accompanying paper by 
iLerchster et al.l (l2010i Paper II hereafter) provides a more de- 
tailed view on the cluster's galaxy population and the performed 
weak lensing analysis. 

The cluster XMMU J1230.3+1339 occupies a special note- 
worthy location in the sky, as it is located behind our closest 
neighbor, the Virgo system. The cluster center is in the direct 
proximity of the Virgo member NGC4477, about 1.3 deg North 
of M87, corresponding to a projected distance of approximately 
400 kpc, or half of R500, at the Virgo redshift. 

This paper is organized as follows: Sect. [2] introduces the 
observations, data reduction, and first basic cluster properties; 
Sect.|3] contains the derived results and mass proxies; Sect.|4] 
focusses on a global pan-chromatic view and the identification 



of cluster associated sub-components; Sect.|5] discusses the dy- 
namical state of the cluster and the comparison of mass prox- 
ies; we conclude in Sect. [6] Throughout the paper we assume 
a ACDM cosmology with Q m = 0.3, = 0.7, and h = 
//o/(100kms _1 Mpc -1 ) = 0.7. Unless otherwise noted, the no- 
tation with subscript X^oq (X200) refers to a physical quantity X 
measured inside a radius, for which the mean total density of 
the cluster is 500 (200) times the critical energy density of the 
Universe p a (z) at the given redshift z. All reported magnitudes 
are given in the AB system. At redshift z = 0.975 the lookback 
time is 7.61 Gyr, one arcsecond angular distance corresponds to 
a projected physical scale of 7.96 kpc, and the cosmic evolution 
factor has a value of E(z) = H(z)/H = [Q m {l + z) 3 + Q A ] 1/2 = 
1.735. 



2. Observations, data reduction, and basic cluster 
properties 

The cluster of galaxies XMMU J1230.3+1339 (see Fig.Q} 
was discovered as part of the XMM-Newton Distant Cluster 
Project (XDCP), a serendipitous archival X-ray survey focussed 
on the identification of very distan t X-ra y luminous system s 
(iMullis et all 120051: iBohringer et ail 120051: iFassbenderl 120071) . 
The first step of the survey strategy is based on the initial de- 
tection of serendipitous extended X-ray sources as cluster can- 
didates in high galactic latitude XMM-Newton archival fields. 
Extended X-ray sources lacking an optical cluster counterpart 
in Digitized Sky Survey (DSS) images are followed-up as dis- 
tant candidates by optical two-band imaging with the objective 
to confirm a galaxy overdensity and obtain a redshift estimate 
based on the color of early-type cluster galaxies. Subsequent 
spectroscopic observations of photometrically identified systems 
provide the final secure cluster identification as bound objects 
and accurate redshifts for the high-z systems. The presented 
XMM-Newton and the VLT/FORS 2 imaging and spectroscopy 
data are part of this XDCP survey data set. 

Additional follow-up data for a more detailed characteri- 
zation of the cluster were obtained with the Large Binocular 
Telescope (LBT) for deep multi-band imaging . Moreover, high- 
resolution X-ray data of a field serendipitously covering the clus- 
ter has recently become available in the Chandra archive. An 
overview of the full multi-wavelength data set used for this work 
is given in Table|2] We first present the optical imaging and spec- 
troscopy data and then discuss the analysis of the archival X-ray 
data. 



2.1. Optical follow-up observations 
2.1.1. VLT Imaging 

Short-exposure images in the R spe ciai (960 s, hereafter R) and 
ZGunn (480 s, hereafter z) filters were acquired in March 2007 
with VLT/FORS 2 in photometric, sub-arcsecond seeing con- 
ditions in order to identify the optical counterpart associated 
with the extended X-ray source XMMUJ1230.3+1339. The 
FORS 2 imaging data with a field-of-view of 6.8'x6.8' were re- 
d uced in a standard m anner following the procedure described 
in lErben et all d2005l) . The 4 (3) bias subtracted and flat-fielded 
frames in the z (R) band were registered to a com mon coordinat e 
system and co-added using SWarp and ScAMP (lBertinll2006l) . 
The final stacked images have a measured on-frame seeing of 
0.52" in the z-band and 0.62" in R. 

Photometric zero-points were derived from dedicated stan- 
dard star observations in the R-band and stellar SDSS pho- 
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Fig. 1. Optical appear- 
ance of the galaxy cluster 
XMMUJ1230.3+1339 
at z = 0.975. Top left: 
Co-added 9.5' x 9.5' i'z' 
image showing the cluster 
environment. The dashed 
circles represent the pro- 
jected cluster radii R200 
(outer) and R 50 o (inner), the 
physical angular scale at 
the cluster redshift is given 
in the upper right corner, 
and the Virgo galaxy 
NGC4477 is marked by 
the dashed circle on the 
right. Top right: RGB color 
composite based on the co- 
added LBT images in i'z' 
(red channel), Vr' (green), 
and UB (blue). The side 
length of the image cor- 
responds to a physical 
scale of 2 Mpc x 2 Mpc 
(250" x 250"), displaying 
approximately the region 
inside Rioo- The large 
number of red cluster 
galaxies are easily visible. 
Bottom left: Core region 
of the cluster in i'z' with 
Chandra X-ray contours 
overlaid in cyan and a side 
length of 560 kpc (70"); 
the green dashed circle 
marks the BCG to the 
South-West of the nominal 
X-ray center (green cross). 
Bottom right: Same sky 
region as on the left as 
color composite. 



Table 1. Observation log of the X-ray, optical imaging, and spectroscopic data coverage of XMMU J1230. 3+1339. 



Observatory 


Instrument 


Data Type 


Date 


Exposure Time 


Observation ID 


XMM-Newton 


EPIC 


X-ray 


Jun 2002 


14ksec/3.9h 


0112552101 


XMM-JVcwton 


EPIC 


X-ray 


Jul 2002 


13ksec/3.6h 


0106060401 


Chandra 


ACIS-S 


X-ray 


Apr 2008 


38ksec/10.6h 


9527 


VLT 


FORS2 


R z imaging 


1 1 Mar 2007 


1.4ksec/0.4h 


078.A-0265 


LBT 


LBC 


UB Vr' i' z' imaging 


28 Feb - 2 Mar 2009 


37ksec/ 10.3 h 




VLT 


FORS2 


MXU spectroscopy 


27 Apr & 6-7 Jun 2008 


2.2ksec/3.9h 


081.A-0312 



tometry in the science field for the z filter. Photometric cat- 
alogs were extracted from the PS F matched images using 
SExtractor dBertin & Arnoutslll996t) in dual image mode with 
the unconvolved z-band frame as detection image. The 50% 
completeness limits^ for the two bands were estimated to be 
z lim (50%) -24.7 mag and R, im (50%) -25.4 mag. 

For this work we limit the quantitative galaxy color assess- 
ment in the cluster environment to this initial FORS 2 R- and 
z-band discovery data set. We constructed the color-magnitude- 
diagram (CMD) for the observed field using the total (MAG_AUT0) 
z-band magnitude and fixed 2.2" (3.5x seeing) aperture magni- 
tudes for the R-z color, both corrected for effects of galactic 



1 The 50% completeness corresponds approximately to a 5 cr detec- 
tion limit (2" apertures). 



extinction. The resulting CMD is displayed in Fig. [2] after the 
removal of stellar sources (CL_STAR>0.95) and objects with sat- 
urated cores. The richly populated cluster red-sequence, i.e. the 
locus of early-type cluster galaxies in the CMD, clearly stands 
out from the background at a color of R-z~2 and provided 
an early indication that the galaxy cluster is indeed at z > 0.9 
prompting the subsequent spectroscopic follow-up. 

We have not attempted to fit a non-zero slope to the clus- 
ter red-sequence, since (i) the slope is expected to be small, (ii) 
the moderately deep ground-based data does not allow precision 
photometry down to faint magnitudes, and (iii) the lack of se- 
cure cluster membership and SED-type information along the 
ridgeline would bias the result due to possible inclusions of in- 
terlopers and non-passive galaxies. For a robust color selection 



4 



R. Fassbenderetal.: A pan-chromatic view of XMMU J1230.3+1339 at z = 0.975 



of predominantly cluster galaxies with a high contrast with re- 
spect to foreground and background objects we apply a simple 
color-cut of +0.2 mag (dotted lines in Fig.|2j around the median 
bright-end color of R-z~2.05 (red dashed line). This color cut 
encompasses 11 out of 13 spectroscopically confirmed cluster 
members and should hence be a good approximation for the un- 
derlying physical red-sequence of the early-type cluster mem- 
bers, which is expected to be slightly bluer at the faint end driven 
mainly by a decreasin g average metallicity with increasing ap- 
parent magnitude (e.g. lKodama & Arimotofl997l) . 

The lower panel of Fig.[2]shows the histogram of the number 
of objects along this color-selected strip 1.85<R-z<2.25mag in 
bins of 0.4 mag for galaxies within 1' of the X-ray center (blue 
solid line) and galaxies further away (black dotted line), with 
the nominal 50 % completeness limit at this color indicated by 
the red vertical line. The apparent decline of the number of red 
ridgeline galaxies towards fainter magnitudes seen in the blue 
histogram can be attributed to several possible factors, or combi- 
nations of these: (i) the increasing uncertainty for fainter galax- 
ies in the R-z color determination scatters a fraction of red clus- 
ter galaxies outside the color cut, (ii) fainter ridgeline galaxies 
are preferentially located beyond a cluster-centric distance of 1' 
(dotted histogram), or (iii) the cluster red-sequence exhibits a 
physical deficit of f aint galaxies as has been reported for seve ral 
high-z clusters (e.g. lDe Lucia et al . 2007l lTanakaetal.ll2007l) . 

Figure[3] displays the cumulative radial distribution of red 
galaxies with the same color cuts as a function of cluster-centric 
distance (blue solid line). About 20 red galaxies are found be- 
yond a cluster-centric distance of 1' but still within the fiducial 
cluster radius (solid vertical line, see Sect. 13.1. Ti l. These galax- 
ies can mostly be attributed to infalling groups, as discussed in 
Sect. l4.2.2l and could contribute to the faint end of the red pop- 
ulation. The blue dotted line shows the average surface density 
of red galaxies within the cluster-centric distance r, which is di- 
rectly related to the galaxy contrast in the enclosed area with re- 
spect to the background. This contrast can be optimized for the 
initial confirmation of the cluster and the identification of its red- 
sequence by restricting the galaxy selection to the inner 30"or 1' 
(vertical dashed lines), as was used for the CMD of Fig. [2] 

For this work, only the particularly richly populated bright 
end of the red ridgeline i n XMMU J1230 3+1339 is highlighted 
(for comparison see e.g. iMei et al 1 120091 A more quantitative 
analysis of the physical cluster red-sequence (and a possible 
truncation thereof) based in the deeper LBT imaging data is dis- 
cussed in Paper II. 

2.1.2. LBT Imaging 

In order to obtain a more detailed view of the cluster, we initiated 
deep wide-field imaging observations in six-bands UB Vr'i' zQ 
with the Large Binocular Cameras (LBC) working in parallel at 
the 2x8.4 m Large Binocular Telescope. The data were acquired 
between 28 February and 2 March 2009 under good, mostly sub- 
arcsecond, seeing conditions. The net exposure times (per single 
8.4 m telescope) over the field-of-view of 26'x26' were 9ksec 
in the U and z'-bands, 5.9ksec in B and i', 3.6ksec in V, and 
3.2 ksec in the r'-band. 



2 For our practical purposes, the difference between the photometric 
reference systems z and z' is negligible (Sj 0.02 mag). The derived LBT 
z' and FORS 2 z-band magnitudes have been confirmed to be consistent 
within the photometric errors. However, for easier distinction we refer 
to z-band magnitudes for the FORS 2 results and z' magnitudes for the 
deeper LBT data. 



The LBT/LBC imaging data were reduced in standard man- 
ner similar to the treatment in the previous section. For more 
information on LBC spe c ific re duction procedures in general 
we refer to iRovilos et alJ d2009l) . and for details on the photo- 
metric calibration and object catalogues based on this data set 
to Paper II. The reduced final co-added image stacks have mea- 
sured seeing full-width half-maxima (FWHM) of 0.8" in the z'- 
band and about 0.9" in all other bands, resulting in limiting mag- 
nitudes (5cr detection in 1" apertures, see Table 1 of Paper II) of 
the co-added frames of: Ui; m ^27.0, Bn m ^26.9, Vii m -26.3, 
r'i; m ^26.2, i'i; m ^25.8, and z'i; m ^25.4. 

This rich imaging data set will be fully exploited in forth- 
coming papers (e.g. Paper II). The current work uses the deep 
wide-field LBT imaging data for two main purposes: (i) a quan- 
titative photometric analysis in the z'-band (Sect. 13.2] ), which is 
deeper than the FORS 2 data and allows accurate background 
determinations in independent external regions; and (ii) a deep 
optical view in single-band and color composites as shown in 
FigsJTJ [lOj Q~TJ and [13] In order to maximize the depth of the 
available data for visualization purposes, we co-added the i' and 
z'-band using an inverse variance weighting scheme to a result- 
ing combined i'z'-frame. This combined band is used as single- 
band background image with a high contrast of cluster galaxies 
and as input for the red-channel of RGB color composite images. 
Accordingly, a co-added Vr'-frame was obtained as input for the 
green-channel, and a UB-image for the blue-channel. 

2.1.3. VLT Spectroscopy 

For the final redshift determination of the cluster, we ob- 
tained spectroscopic observations with VLT/FORS 2 using a sin- 
gle MXU-mode (Mask eXchange Unit) slit-mask centered on 
XMMU J1230.3+1339 for a total net exposure time of 7.85ksec 
(TableQ}. Targeted galaxies were selected to be preferentially 
red galaxies close to the ridgeline in Fig|2] The chosen instru- 
ment setup with the 3001 grism and a slit width of 1" provides 
a wavelength coverage of 6 000- 10 000 A with a resolution of 
R = 660. The data consisted of six individual exposures taken in 
three different nights under variable seeing conditions ranging 
from 0.7" to 1.9". 

The spectra were reduced with standard techniques using 
tasks. In short, the 2-D spectra were bias subtracted, flat- 
fielded to correct for pixel-to-pixel variations, and co-added to 
provide a sufficient signal-to-noise ratio (SNR) for the tracing of 
weak spectra. 38 1-D spectra were extracted and sky-subtracted 
by fitting a low order polynomial to adjacent background regions 
along the trace of the spectrum. The wavelength calibration was 
achieved by means of a Helium- Argon reference line spectrum 
observed through the same MXU mask. The quality of the ob- 
tained dispersion solutions was cross-checked by measuring the 
observed position of several prominent sky lines yielding typical 
absolute rms calibration errors below 1 A. 

For the redshift determination, the reduced and wavelength 
calibrated spectra were cross-correlated with a r ange of galaxy 
templ ate spectra using the IRAF package RVSAO dKurtz & Minkl 
1 1998b . Out of 38 spectra, we could identify 13 secure clus- 
ter members (see Tabled with a median cluster redshift of 
z = 0.9745 ± 0.0020 and a bootstrapped error estimate. Figure|4] 
displays five sample spectra of cluster members with the most 
prominent spectral features labelled. The cross-correlation of 
absorption-line spectra with a sufficient signal-to-noise ratio 
(SNR^5) yields typical statistical velocity errors of ±50km/s 

3 http://iraf.noao.edu 
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Fig. 2. Top: R-z color-magnitude diagram of the cluster envi- 
ronment based on the FORS 2 imaging data. The richly pop- 
ulated cluster red-sequence around the median ridgeline color 
R-z~2.05 (red dashed line) is clearly visible. Filled red cir- 
cles indicate objects within a 30" radius from the X-ray center, 
green circles correspond to 0.5-1', small black dots represent 
all other objects in the field, and open circles mark spectroscop- 
ically confirmed cluster members. Horizontal dotted lines de- 
fine the applied color cut interval of +0.2 mag about the median 
color, while the dotted vertical line refers to the expected appar- 
ent magnitude (m*^2 1 .9 mag) of a passively evolving L* galaxy 
at the cluster redshift (in both panels). The 50% completeness 
limit is indicated by the black solid line. Bottom: Histogram of 
galaxy counts in the red color interval 1 .85<R-z<2.25 in bins of 
0.4 mag for galaxies within 1' from the cluster center (blue line). 
The vertical solid line indicates the nominal 50% completeness 
limit of the data, while the black dotted line shows the magnitude 
distribution for galaxies at larger distances for comparison. 



for c-z. Adding the wavelength calibration uncertainty (£1 A at 
8 000 A) in quadrature, we estimate the typical uncertainty for 
most of the secure individual galaxy redshifts in Table|2]to be 
o- z ~ 0.0002 

Due to the partially poor seeing conditions, not all redshifts 
of targeted galaxies could be determined with high confidence. 
We have identified 1 1 additional tentative galaxy redshifts which 
are close to the median cluster redshift and potentially indicate 
member galaxies. For completeness these objects are also dis- 
played in Fig. [5] which shows the restframe velocity distribu- 
tion as a function of cluster-centric distance further discussed in 
Sect. l3~2~6l 
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Fig. 3. Cumulative radial distribution of color selected red galax- 
ies (as in Fig.[2j as a function of cluster-centric distance (blue 
solid line). Dashed vertical lines indicate the radii used for color 
coded symbols in Fig. [2] the solid line marks the fiducial clus- 
ter radius. The dotted blue line shows the cumulative radial 
red galaxy number counts normalized to the enclosed area, i.e. 
the average red galaxy density per square arcminute within the 
cluster-centric distance r. 



1 .0 
0.5 
0.0 

0.5 

0.0 

1.0 

0.5 
0.0 

1.0 

0.5 
0.0 

1.0 

0.5 
0.0 

6000 6500 7000 7500 8000 
observed wavelength (A) 



z = 0.975 


j ] 

I , : 


( \ \ i 

h L II I VV 

/\/N i ii 

.i.i. . i . . . ii. 




Oil 

J ■ 

A/^^W | Vlj 


K H H6 

Ii ifU^VWvV 

Mv! 1 
J 11 W 1 

1 .1 


, GHy 

VW IV 
1 l 
l 




'■ 
'■ 

<J*\j^ |^^ 


1 l k l l fn 

i.uAvu /VyVi/Vy 

AtAr r*^ r i - 

/ V i i 

1 M 




'■ 


ii i i i jA. 

fififl 1 1 

vy i ii 

— U — ^-1 1-4— 




■ 


1 1 1 lij 
i (|AmIn^^Ii/j^ r 

a'a'J 1 I 

Vv , 1 (I .. II 



8500 



Fig. 4. Sample spectra of secure cluster members with the me- 
dian redshift of z — 0.975, smoothed with a 9 pixel boxcar filter. 
The expected observed positions of prominent spectral features 
at the median redshift are indicated by blue dashed lines. The 
sky spectrum (bottom) and telluric features (top) are overplotted 
in red. From top to bottom the spectra correspond to the galaxies 
with object IDs 1-9-4-6-8 in Table|2] 



2.2. X-ray Data 
2.2.1. XMM-Newton 

XMMU J1230.3+1339 was initially detected as an extended X- 
ray source at high significance level in the XMM-Newton field 
with observation identification number (OBSID) 0112552101 
at an off-axis angle of 4.3' from the central target, the nearby 
galaxy NGC4477 in the Virgo cluster. The original XDCP 



6 



R. Fassbenderetal.: A pan-chromatic view of XMMU J1230.3+1339 at z = 0.975 



3000 



2000 



1000 



-1000 



-2000 



-3000 



□ • 



□ 







200 



400 



600 800 
R [kpc] 



1000 1200 



1400 



Fig. 5. Cluster restframe velocity distribution of spectroscopic 
cluster members as a function of projected distance from the X- 
ray center. Member galaxies of Table|2] with secure redshifts are 
indicated by blue circles, dashed lines correspond to the formal 
velocity dispersion estimate. For completeness, galaxies with 
tentative redshifts are also shown as squares, smaller black sym- 
bols indicate velocities outside the formal 3 <x interval. 

Table 2. Secure spectroscopic cluster members of 
XMMU J1230.3+1339 with total z'-band magnitudes, cluster- 
centric distances, and measured redshifts. 
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Fig. 6. Growth curve analysis for the unabsorbed, background 
subtracted 0.5-2 keV X-ray flux as a function of projected cluster 
radius. The blue solid line shows the total 0.5-2 keV flux inside 
the projected radius r as measured with the XMM-Newton PN 
camera, the red line represents the sum of the two MOS instru- 
ments. Poisson errors plus 5% uncertainties in the background 
determination are displayed by the dashed lines, the vertical line 
depicts the R500 radius. The fully consistent Chandra GCA re- 
sult is overplotted in black. The X-ray emission can be traced 
out to about 90", where it reaches the plateau value. 



ID 


RA 


DEC 


z' 




£spec 




J2000 


J2000 


mag 


arcsec 




01 


12:30:16.41 


+ 13:39:16.2 


20.37 


15.1 


0.9720 


02 


12:30:16.35 


+ 13:38:49.5 


19.90 


18.0 


0.9787 


03 


12:30:18.43 


+ 13:39:06.6 


21.00 


21.1 


0.9794 


04 


12:30:17.76 


+13:39:26.0 


20.98 


24.4 


0.9758 


05 


12:30:20.11 


+ 13:39:00.9 


21.89 


46.3 


0.9690 


06 


12:30:20.44 


+ 13:38:38.6 


20.97 


57.2 


0.9767 


07 


12:30:14.95 


+ 13:39:58.3 


21.57 


62.3 


0.9763 


08 


12:30:18.87 


+ 13:38:00.0 


20.09 


69.9 


0.9742 


09 


12:30:11.94 


+13:40:33.4 


20.96 


117.3 


0.9692 


10 


12:30:09.58 


+13:40:38.4 


23.00 


146.1 


0.9674 


1 1 


12:30:21.50 


+ 13:36:54.1 


21.86 


146.4 


0.9745 


12 


12:30:09.16 


+13:40:49.3 


21.40 


158.0 


0.9658 


13 


12:30:06.11 


+ 13:40:22.7 


22.81 


181.6 


0.9745 



source detection run, for which 470 XMM-Newton archival 
fields were processed, was performed with SAS v6.5 using the 
tasks eboxdetect and emldetect. With an extent likelihood of 
Lext -117, corresponding to a formal probability that the source 
extent is due to a statistical Poissonian background fluctuation 
of ppoisson - 1.5 x 10" 51 , the cluster is one of the most signifi- 
cant X-ray sources among all distant cluster candidates in the 
XDCP survey. A detailed descriptio n of the source dete ction and 
screening procedure can be found in lFassbendej d2007h 

For the analysis presented here we re-processed the field 
with the latest version SASv8, and additionally considered a 
second archival field with OBSID 0106060401 (Tabled), which 
encloses the cluster source at a larger off-axis angle of 9'. 
After applying a strict two-step flare cleaning process for the 
removal of high background periods, a clean net exposure 
time of 13.3/8.3 ksec remained for the EMOS/PN cameras in 
field 0112552101, and 7.8/5.0ksec for the second observation 
0106060401. 

In order to optimize the signal-to-noise ratio of low surface 
brightness regions in the cluster outskirts discussed in Sect. 14.21 



we combined the X-ray images of both XMM-Newton observa- 
tions and their corresponding exposure maps with the SAS task 
emosaic. This results in more than 1 100 source photons from 
XMMU J1230.3+1339 in the mosaic image for the used 0.35- 
2.4 keV detection band , which provi des an optimized SNR for 
distant cluster sources (IScharfll2002h . To enable the evaluation 
of the significance of low surface brightness features in the pres- 
ence of chip gaps, we created a signal-to-noise map of the clus- 
ter environment. Using this map, we created a log-spaced con- 
tour plot spanning the significance levels 0.5<SNR<35 used in 
Fig.[[0]of Sect.gS 

We applied the g rowth curve analysis (GCA) method of 
iBohringer et al.l d2000l) to obtain an accurate 0.5-2 keV flux mea- 
surement of the cluster as a function of radius. Since the effec- 
tive exposure time and the point-spread-function (PSF) at the 
cluster source are significantly better in field 0112552101, we 
focus the quantitative analysis on this single observation and 
use the second field 0106060401 only for a consistency check. 
Contaminating regions containing X-ray sources not associated 
with the cluster were conservatively masked out prior to the 
GCA measurement. The X-ray emission of the neighboring SO 
galaxy NGC4477 (mostly point sources) only contributes signif- 
icantly at cluster-centric radii beyond 2', i.e. the X-ray analysis 
of XMMU J 1230.3 +1339 is basically not influenced by the fore- 
ground galaxy. 

Figure|6] displays the GCA results for the cluster emission, 
which can be traced out to a radius of about 90", i.e. just beyond 
the nominal radius R500 - 84" as estimated in Sect. 13. II Within 
the projected R500 radius we measure an unabsorbed 0.5-2 keV 
flux of f x ,500 = (5.14 + 0.54) x 10~ 14 erg s^'crrT 2 . For complete- 
ness we also show the fully consistent Chandra GCA result in 
Fig.|6](black line). 

For the temperature determination of the cluster we f ollowed 
the general procedure described e.g. in lPratt et al.l d2010t) . As for 
the flux determination, we focused the spectral analysis on field 
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Channel energy (keV) 

Fig. 7. Chandra and XMM-Newton X-ray spectra of 
XMMUJ1230.3+1339 with a joint fit ICM temperature of 
7x,2500 - 5.3 keV. The Chandra spectrum extracted from the 
ACIS-S3 chip is shown in blue, and the corresponding data 
from the three XMM-Newton instruments in green (PN), black 
(MOS1), and red (MOS2). 



0112552101, since the effects of the broadened XMM-Newton 
PSF at larger off-axis angle outweigh the gain of additional 
source photons. We extracted an X-ray spectrum of the source 
from two circular regions of radius 40" and 71" centered on 
the X-ray peak. The smaller 40" aperture is signal-to-noise op- 
timized in conjunction with the Chandra spectral analysis, the 
larger aperture corresponds to the maximum radius for which re- 
liable results can be derived with XMM-Newton^ The presence 
of extended thermal foreground emission from the Virgo cluster 
results in an increased background leveQ, which we assume to 
be homogeneous over the cluster scale of interest of about 2'. 
To account for this extra foreground component, we have tested 
several methods with both local and external background spec- 
trum determinations. The most robust results were obtained for 
the 40" aperture using a local background spectrum extracted 
from a nearby uncontaminated region after the masking of point 
sources. With the metal abundance fixed to Z=0.3Z o , we derive 
a best fitting core region temperature for XMMU J1230.3+1339 
of 5.28+°;^ keV, based on a single temperature MEKAL model, 
a minimally binned spectrum (>1 cts/bin), and C-statistics. The 
XMM-Newton spectral fit and the residuals are displayed in 
Fig. [7] (green, red, and black lines). 

For the larger 71" aperture, we measured a slightly higher 
X-ray temperature of 5.6+j A , keV applying the same local back- 
ground subtraction procedure. Since this larger extraction region 
exhibits a lower signal-to-noise ratio and is hence more prone 
to background uncertainties, we cross-checked this temperature 
trend with an alternative spectral fitting method using an ex- 

4 The source-flux-to-total-background-fiux ratio (0.5-2 keV) is about 
3:2 for the 40" aperture and 1:2 for the 71" aperture. 

5 The Virgo contribution to the total background is about 80%. 



ternal background spectrum. In addition to the physically mo- 
tivated model of the cosmic X-ray background (CXB) consist- 
ing of two unabsorbed MEKAL models p l us a power law with 
fixed index y = 1.4 (see iLumb et al1l2002t iDe Luca & Molendil 
l200l . we considered a further thermal MEKAL component 
of the total background spectrum in order to account for 
the local foreground emission of Virgo at the position of 
XMMU J 1230. 3+ 1339. From this four-component model we 
obtained a best fitting local Virgo temperature of kTvi rgo = 
1.88 ke V consistent with p revious ASCA measurements of this 
region dShibata et al.ll200lt) . This four-component model, with a 
renormalization appropriate to the ratio of the surface area of the 
extraction regions, was then added to the fit of the cluster spec- 
trum, resulting in an X-ray temperature based on this external 
background model of 6.4+] 2 keV. 

2.2.2. Chandra 

Additional X-ray coverage of XMMU J1230.3+1339 has re- 
cently become publicly available in the Chandra archive, where 
the cluster was serendipitously observed almost on-axis in a 
38.2 ksec pointing targeted on NGC4477 (Observation ID 9527, 
see Tab.Q]). The field was observed in VFAINT mode, with the 
cluster being located on the central ACIS-S3 detector at an off- 
axis angle of approximately 2'. 

We performed a standard data reduction using the CIAO 
v4 . 1 software package, with a recent version of the Calibration 
Database (CALDB v4 . 1 . 1). Since the observation was taken in 
VFAINT mode, we ran the task acis_process_events which 
provides an improved flagging of probable background events by 
distinguishing real X-ray photon events and events most likely 
associated with cosmic rays. With this procedure the ACIS par- 
ticle background can be significantly reduced. The data is fil- 
tered to include only the standard event grades 0, 2, 3, 4 and 
6. We confirmed that there were no flickering pixels with more 
than two events contiguous in time and checked visually for hot 
columns remaining after the standard cleaning. The resulting to- 
tal net exposure time after this cleaning procedure is 37.68 ksec. 

We repeated the growth curve analysis for Chandra in the 
soft 0.5-2.0 keV band (Pig.|6j, which traces the cluster emission 
out to similar radii as in the single XMM-Newton field. After 
subtraction of the background of b- 0.043 +0.008 cts/pixel mea- 
sured in an external region 634 net cluster counts within r=2' re- 
main. The resulting unabsorbed Chandra fluxes are in full agree- 
ment with the specified XMM-Newton values as shown by the 
black line in Fig. [6] 

The main advantage of the high-resolution Chandra X-ray 
data is the ability to resolve small-scale ICM structures in the 
inner cluster region and to obtain precise radial surface bright- 
ness (SB) profiles down to the central core. We measured the 
binned azimuthally-averaged surface brightness profile and per- 
formed a fit with th e widely used approximation of th e single 
isothermal ft model (ICavaliere & Fusco-Femianoll 19761) . tracing 
the ICM out to about R500. We fitted the functional form S (r) = 
S (l + ir/r c ) 2 T 3l3+0 - 5 + b for the radial SB profile Sir) with free 
parameters So for the central SB, r c for the core radius, ft for the 
slope, and a constant b for the background using a Levenberg- 
Marquardt least-squares minimization. The radial profile and 
the best fitting single ft model are shown in the left panel of 
Fig. [8] with resulting parameters So — 0.28 + 0.04cts/pi:xQ, ft = 
0.837 + 0.471, r c = 215 + HOkpc, and a reduced-chi square of 
0.70. Beyond a projected radius of about 150 kpc, the model is a 



Using 2x2 binned pixels with an effective scale of 0.984" per pixel. 
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Fig. 8. Chandra surface brightness profile of XMMU J1230. 3+1339 (black data points) with the best fitting single (left) and two- 
component (right) beta models (red solid lines). The background level is indicated by the blue line. The dotted lines in the right 
panel represent the individual fit components, the inner one includes the background and the outer one is additive. 



very good representation of the data. The reasons for the slight 
deviations in the core region at r 55 50 kpc and the fiat shoulder 
at r ~ 100 kpc can be understood from a more detailed look on 
the cluster core region (Sect 13. and the dynamical structure 
in Sect.E3]and Fig. [TBI 

In order to obtain an improved model representation for the 
innermost bins of the radial profile we also performed a two- 
component p model fit, with a single /3 coefficient, but allowing 
for two independent core radii r c \ and r c i 



S(r)=Si 



1 + 1 — 



-3jS+i 



+ s 2 



1+1^ 

r C 2 



-3j8+i 



+ b 



The resulting fit parameters for this model are f3= 1 .0 (the maxi- 
mum allowed value), Si =0.42 ± 0.51 cts/pix, r c \ = 42 ±34 kpc, 
S 2 = 0.24 + 0.04 cts/pix, and r c i -212 + 32 kpc with a reduced- 
chi square equal to 0.45 As shown in the right panel of Fig.[8j 
the core region of the SB distribution at r $ 100 kpc is now qual- 
itatively better modelled. 

We also measured the surface brightness concentration for 
the cluster using the para meter csb = ( S -Bfr < 40kpc)])/(SB[r < 
400 kpc]) as defined by ISantos et al.1 d2008h . which allows an 
improved structural characterization of high-z clusters com- 
pared to methods requiring very high photon statistics. For 
XMMU J1230.3+1339 this parameter is found to be c SB = 0.07, 
which places the system in the formal category of Non-Cool 
Core (NCC) clusters, in contrast to the more centrally concen- 
trated objects in the moder ate and strong Co ol Core Cluster 
(CCC) class at high redshift dSantos et al.ll2008l) . 

For the visual representation of the two dimensional X-ray 
surface brightness structure (Figs. 11191 and [T3l , we adaptively 
smoothed the Chandra soft band image and created an X-ray 
contour map with eight log spaced surface brightness levels 
spanning the significance range of 1 -70 cr above the background 
rms. It should be noted that the cluster does not contain any de- 
tectable point sources within Z?5oo at the current Chandra depth. 
The structural and spectral analysis of the cluster properties 
should thus be unbiased, also for XMM-Newton with its lower 
spatial resolution. 



We performed a spectral temperature fit in a signal-to-noise 
optimized apertur e of 40" following the procedure described in 
iTozzi et all (120031) . Applying the same general procedure as for 
XMM-Newton with the local background measured in an on- 
chip region with 60" diameter, we obtain a best fitting Chandra 
core temperature of 5.40^ ^keV with the abundance fixed at 
0.3 Z & . The spectrum and model fit is displayed in blue in Fig. [7] 
which also shows an indication of the Fe line at kT ~ 3.3 keV. 
Leaving the metal abundance as a free parameter, we can derive 
a weak Chandra-only core region constraint of Z = 0.46Q2g 8 Z . 



(!) 3. Results and mass estimates 



In the following section we will present results derived from the 
X-ray measurements and the optical imaging and spectroscopic 
data. A special focus will be given to a range of observables 
that are commonly used as proxies for the cluster mass. All final 
M200 cluster mass estimates (ME) based on a given observable 
are labelled with a number (e.g. ME 1) for the purpose of easier 
referencing, intermediate results are indicated by an additional 
letter (e.g. ME la, ME lb). All used mass-observable scaling re- 
lations are stated explicitly for clearness, transformed into our 
assumed cosmology and the used cluster mass definition M200, 
where applicable. 

3. 1 . Global X-ray properties 

We will start with the X-ray determined bulk properties of 
the intracluster medium in order to establish some first ro- 
bust fundamental characteristics of the cluster. For all subse- 
quent discussions, we define the 'center-of-mass' of the ICM X- 
ray emission determined with XMM-Newton (RA=12:30:16.9, 
DEC=+13:39:04.3, see Tab.0 as the reference center of the 
cluster, and measure all cluster-centric distances relative to this 
point. 



3.1 .1 . X-ray temperature 

The Chandra and XMM-Newton results on the ICM tempera- 
ture in the 40" aperture core region are of comparable statis- 
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tical quality and are fully consistent with each other. In order 
to further decrease the uncertainties, we have performed a joint 
Chandra /XMM-Newton spectral fit to these data for the region 
within 40" = /?2500 (see below). We find a final joint X-ray tem- 
perature of kT x4Q // ^ kT x ,2500 = 5-30^ g2keV (see Fig.[7j us- 



O.43 + °;^Z , i.e. with a 2.3 cr 



ing the same method as in Sect. 12. 21 with a fixed metallicity of 
Z=O.3Z , minimally binned data, and C-statistics. Leaving the 
metallicity as a free parameter, we can constrain the iron abun 
dance in the cluster core to Z2500 
detection significance. 

The XMM-Newton data analysis indicates that the global 
ICM temperature of XMMU J1230.3+1339 is higher when the 
measurement aperture is increased. With the currently available 
X-ray data set this is a tentative trend, as the measurement er- 
rors also increase. However, to allow for an evaluation of pos- 
sible biases, we take the average of the 71" temperatures (lo- 
cal and external background) as a second reference temperature 
kT X7[ // =* 6.0^2 keV. Since ideally one would like to have a 
measurement of Tx,5oo, i.e. within the full R500 region, the latter 
value in the 6 ke V range could indeed be a better approximation 
of the underlying global temperature of the cluster. 

This global X-ray temperature is the fundamental parame- 
ter of the intracluster medium and is closely linked to the depth 
of the underlying potential well of the cluster through the Virial 
theorem. We use the jointly determined best (core) X-ray tem- 
perature Tx,2500 with an uncertainty of about +13% to derive 
robust size measures of the cluster, which is an important pre- 
requisite for the measurement of most clus ter related quantities . 
We are using the R-T scaling relations of lArnaud et al.l (120051) . 
which have been calibrated using a local cluster sample based on 
high-quality XMM-Newton dataj. Using the R-T scaling relation 
for hot systems (r x > 3.5 keV) R 2Q0 = (1714 ± 30) x E~ l (z) x 
(kT/5 keV) a50±0 05 kpc, we obtain for XMMU J1230.3+1339 an 
estimated /?2oo = 1017^ kpc. For R500, using the normaliza- 
tion 1129 + 17 kpc of the R-T relation, the estimated radius is 
#500 = 670^ kpc, and similarly 7?25oo = 297+ j 7 kpc (normal- 
ization 500 + 5 kpc). Standard error propagation was applied in 
all cases to evaluate the uncertainties. In terms of projected dis- 
tance in the plane of the sky at the cluster redshift, these repre- 
sentative length scales correspond to angular distances of [R200, 
#500, fl 2 50o]=[127.8", 84.2", 37.3"]. With respect to the higher 
temperature estimate T x 71 », all derived radii would increase by 
about 6%, i.e. possible systematic effects are mostly small in 
comparison to other uncertainties. 

As the next step, we can derive a first mass es t imate of the 
cluster by using the M-T relation of lArnaud et alJ d2005h . cali- 
brated with the same local cluster sample as above. Using again 
the best joint temperature Tx,2500, we obtain from the relation 



Tx (5.74 ± 0.30) 

''"200 



E(z) 



kT 



5kev) 



1.49±0.17 



x 10 14 M Q 



(2) 



an estimate of the approximate total mass of the cluster of 
M 200 = 3.61±°^ X 1O 14 M (ME la). For M 500 the M-T rela- 
tion, with the normalization of (4.10 + 0.19) x 1O 14 M0, yields 
M 50 o = 2.58 + ^ 4 xl0 14 M . 

The use of the higher temperature estimate T x 71 « translates 

into a 20% higher mass of M 200 - 4.34 x 10 14 M o (ME lb). 
We take this positive offset of 0.73 x 10 14 Mq into account 



7 Temperatures in lArnaud et alJ d2005f ) are determined in the region 
[0. l-0.5]i?200- The effect of using T X2 5oo f° r this work should be small 
and is accounted for by the discussed systematics. 



as a potential systematic error (denoted in brackets) and arrive 
at a final X-ray temperature-based mass estimate of M^L - 

3.6l!°^ (+O ' 73) xlO 14 M (MEl). 

One of the main objectives of this work is the compari- 
son of total cluster mass estimates, most commonly approxi- 
mated by M200, as derived from various mass-observable rela- 
tions. For scaling relations calibrated to yield only M500, we 
apply the ratio of the normalizations of the respective M-T re- 
lations (see above) /5oo^200 = 200 /A/500 - 1-40 to obtain a 
total mass estimate M2oo- In the same way, /25oo^200 is given by 

/2500^200 = ^200/^2500 - 3.21. 



3.1.2. Hydrostatic mass estimate 

Under the assumption of an ideal isothermal gas in hydrostatic 
pressure equilibrium (HE), the following widely us ed approxi- 
mation for the total mass profile can be derived (e.g. lHicks et al.l 
l2^0lllRosati^t7Dl2009l) 



Of 



G fi m p 



(3) 



1.13 



Mpc 



where ks is the Boltzmann constant, G the gravitational con- 
stant, m p the proton mass, and /i =* 0.59 the mean molecular 
weight in atomic mass units. By using the joint X-ray tem- 
perature T Xj 25oo> an d the single beta model fit parameters from 
Chandra (J3 -0.84, r c ^215 kpc), we obtain the hydrostatic mass 
estimates for XMMU J1230.3+1339 of M 500 ^3.1xl0 14 M o and 
A/200 ^4.9xl0 14 M Q (ME 2a), whereas the latter value is slightly 
extrapolated beyond the last reliable data point of the measured 
X-ray SB profile. 

Note that only about 11% (18%) of the total mass M 20 o 
(M500) is actually enclosed within one core radius r c , within 
which the single beta model exhibits some deviations from 
the data possibly due to the onset of cooling (see Sect. 13.1. 41 
and significant dynamical activity (see Sect. l4.4b . Between one 
core radius and r ^ R500 the model describes the radial sur- 
face brightness profile with good accuracy, so that an extrapola- 
tion out to R200 is justified. The dominant source of uncertainty 
here is the ~56% error on the ft parameter due to a degener- 
acy with the fitted core radius over the accessible data range. 
Taking this and the temperature and radii errors into account 
we arrive at a final hydrostatic equilibrium mass estimate of 
M™ = (4.9 ± 2.8) x 1O 14 M (ME 2). 



3.1.3. X-ray luminosity 

From an observational point of view, the X-ray luminosity Lx is 
the ICM bulk property which is easiest to determine, as it only 
requires flux and redshift measurements. For this reason, the X- 
ray luminosity is of crucial importance for most large-scale X- 
ray surveys and provides t he link to other ICM proper ties and 
the total cluster mass (e.g. iReiprich & Bohringerll2002l) via Lx 
scaling relations. 

For XMMUJ1230.3+1339 we determine a restframe 0.5- 
2 keV cluster luminosity of i^soo^ = (l-92±0.20)x 10 44 erg/s 
within the accessible aperture Rsoo- The bolometric X-ray lumi- 
nosity can be determined by applying a temperature dependent 
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bolometric correction factor f£ Y-2tev to tne so ^ band luminosity 
(e.g. lBohringer et al. 2004]), in our case f£f_ 2 ^ ey -3.38, yielding 
a total bolometric [0. 01-100 keV] luminosity for the cluster of 
^500 = /o.t2 ke v-4,50o keV = (6.50 + 0.68) x 10 44 erg/s. 

Recently the use of a core-excised X-ray luminosity has been 
shown to re duce the scatter in the Lx scaling relations by a 
factor of 2 dPratt et al.l l2Q09h . simply by avoiding the central 
cluster region and its related complex physical processes. For 
XMMU J1230. 3+1339 we measure a core-excised luminosity 
in an aperture of [0.15-l]# 5 oo of ^,[0.i5-T]50O = ( L53 ± °- 21 ) x 

10 4 4 erg/s. 

iPratt et al.l (120091) also provide the latest Lx scaling rela- 
tions calibrate^_wjfhdierepresentative local REXCESS sam- 
ple dBohringer et al.l 120071) . Applying the self-similar L-T re- 



lation E~ l (z) 



rbol 

^X,500 



(7.13 + 1.03)(/tr/5keV) 



3.35±0.32 



10 44 erg/s, we would predict an X-ray temperature for the clus- 
ter of ^4.1 keV, i.e. a value approximately 2cr below our spec- 
troscopic temperature measurement. The main reason for this 
discrepancy is the assumption of self-similar evolution of the 
X-ray luminosity, i.e. Lx oc E{z) at fixed temperature. While 
some empirical result s at z $ 1 are consist e nt with this evolu- 
tionary scenario (e.g. IVikhlinin et al ] l2002t iKotov & Vikhlininl 
12005b iMaughan et al.l l2006h . other studies exten ding towards 
higher redshifts find lit tle or no evolution (e.g. lO'Hara et al.l 
120071 iRosati et al.ll2002l) . or even negative trends with redshift 
(lEttori et all l2004 . XMMU J1230.3+1339 by itself would be 
consistent with either of the latter two scenarios. Since the evolu- 
tionary trends of the X-ray luminosity at high redshifts (z > 0.8) 
are not well constrained yet, we assume as a starting point a no- 
evolution scenario, which would predict Tx - 4.9 keV for the 
given luminosity, and treat a possible Lx redshift evolution as 
a systematic uncertainty. The use of a no-evolution scenario for 
the L-T relation instead of the self-similar prediction is also em- 
pirically supported by the latest compilation of high-z cluster 
measurements (Reichert et al., in prep.). 

With this assumption on the evolution of Lx, we can derive 
an additi onal mass estima te of the cluster based on the L-M re- 
lation of lPratt et al.l (120091) 



rbol 

^X,500 

Ekz) 



(1.38 + 0.12) 



00 



2x 10 14 M 



O 



2.08 ±0.13 



xlO 44 erg/s, (4) 



3.1.4. Core properties of the ICM 

We now have a more quantitative look at the properties of the 
intracluster medium in the densest central core region. The line- 
of-sight emission measure EM = J n 2 dl-^n(\+zf S (r)/A(T, z) 
gives direct observational access to the local electron density « e , 
deduced from the deprojected measured surface brightness pro- 
file S(r) as a function of cluster-centric distance r in the soft 
X-ray band, and the redshifted and absorption corrected cool- 
ing function A{T, z). For an assessment of the central density 
of the cluster, we use the double beta model surface brightness 
fit (right panel of Fig. [8), since the inner beta model component 
dominates the densities within the central 40kpc. After the de- 
projection of the SB model fits (Si, r c i, S2, r C 2, ft) of Sect. l2.2.2l 
we obtain a total electron density of n e o - (2 . 3 + 1 .4) x 1 0~ 2 cm -3 
in the very center, and n e ,20kpc - (1.7 + 1.0) x 10~ 2 cirT 3 at a 
cluster-centric radius of 20kpc. 

The central density and temperature of the ICM determin e 
the cooling time scale f coo i in the cluster core (ISarazi 
i.e. the time for which the thermal energy of the gas would be 
radiated away at the current cooling rate 



f coo i-2.9xl0 lu yr-| 



10- 



\~UkT \ m 
'] \11 



keV/ 



(5) 



From this expression we derive a central cooling time scalJ^l 
for XMMU J1230.3+ 1339 of f coolj0 = (2.9 + 1.8) Gyr and at the 
20kpc radius f CO oi,20kpc = (4-0 + 2.4) Gyr. This time scale is to be 
compared to the age of the Universe at the redshift of the cluster 
(f - 5.85 Gyr), or for more realistic cluster formation scenarios, 
the time span since z = 3 (t 3.74 Gyr) or z = 2 (f =a 2.63 Gyr). 

The central electron density « e ,o provides an alternative 
test for th e presence of a co ol core. Following the CCC def- 
inition of IPratt et al] d2009l) (n e ,o > 4 x 10~ 2 • £ 2 (z)cm 3 ), 
XMMU J1230.3+1339 is consistently classified as a Non-Cool- 
Core cluster, in concordance with the surface brightness concen- 
tration measure csb in Sect. l2.2.2l However, the cluster's core 
region is in an interesting transition phase, where the dense gas 
has basically been going through a single full cooling time scale 
since redshift 2-3. Hence, the central core region could indeed 
contain gas which has already cooled considerably, consistent 
with our tentative finding of an increasing average X-ray tem- 
perature with radius in Sect. l3.1.T1 



corrected for one factor of E(z) to take out the self-similar 
redshift scaling (in L-T). We thus obtain a mass estimate of 
M 500 - (2.96 + 0.41) x 1O 14 M from the L-M relation, which 
yields M 200 ^ / 5 oo^200 • M 500 * (4.1 ± O.6)xlO 14 M (ME 3a). 
However, the dominant uncertainty in this luminosity-based 
mass proxy is the discussed weak constraint on the luminos- 
ity evolution. We take this current lack of knowledge into ac- 
count by adding a 23% systematic error, based on the decreased 
mass estimate of the self-similar evolution model. The final 
luminosity -based total mass estimate is hence — (4.1 ± 

0.6 (+1.0))xl0 14 M o (ME 3). Using the core excised X-ray lumi- 
nosity ^xio^5-T]500 W ^ m tne corres P onc ling scaling relation yields 
a consistent mass estimate within 8%, implying that the central 
region of XMMU J1230.3+1339 does not exhibit any particular 
feature that would drive a deviation. 



8 Derived with the BCES orthogonal fit method and corrected for 
Malmquist bias. 

9 This approach is somewhat ad hoc at this point, but provides a sig- 
nificantly improved agreement with other high-z cluster measurements. 



3.1.5. Gas mass 

The total gas mass M gas of clusters is obtained from the in- 
tegrated radial gas density profile. For XMMU J1230.3+1339, 



,-3 



both the deprojected single beta model (« e ,o - 0.71 x 10 cm 
ft 0.84, r c 215 kpc) and the double beta model yield consis- 
tent results for the gas mass of M gaSi50 o - (3.0 ± 0.9) x 10 13 M@, 
where the estimated error is based on the uncertainties of the 
latter model in combination with the allowed temperature range. 

Using the mean local gas mass fract ion for similarly mas- 
sive clusters of / gas , 5 oo - 0.10 ± 0.02 (IPratt et alJ 120091) . we 
readily obtain a fourth total mass estimate of /5oo^200 • 

M gas , 500 //ga S ,500 - (4.2 + 1.5) x 1O 14 M (ME 4). 

With the knowledge of the final combined best estimate for 
the total cluster mass as discussed in Sect. 15.21 we can 

turn the argument around and estimate the cluster gas mass frac- 
tion for XMMU J1230.3+1339 for the assumed cosmology. This 



10 When considering only the single /3-model fit parameters, this cool- 
ing time scale increases by a factor of 3. 
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consistency check yields / gas , 5 oo - Mga^oo/CM^V/soo-aoo) 
0.100 + 0.035 in good agreement with the average local value. 

3.1.6. Y x 

As the final Chandra and XMM-Newton derived property, we 
determine the X-ray version of the Compton-Y parameter Yx = 
M g - d s ■ Tx- This quantity reflects the total thermal energy of the 
ICM and is predicted to be a robust, low s catter mass proxy 
from simulations (e.g. iKravtsov et al.ll2006t iNagai et aill2007l) 
which has re cently been well confirmed by observations (e. g. 
lArnaud et al.ll2007t l Vikhli nin et al . 20093 lArnaud et al.ll2010l) . 

From the measurements of the global X-ray temperature 
and the gas mass within the radius R500, we obtain lx,500 = 
(1.6 ± 0.5 (+0.2)) X lO 14 M0keV, where again a systematic er- 
ror is considered for a possible temperature bias CSect. 13. lTTb . 
With the latest cali bration of the M-Y^ scaling relation from 
lArnaud et ail (l2010l) we find 

inl4.567±0.010 / y \0.561±0.018 

M500 = - 2 - . * - M (6) 

Et( z ) \2 x 10 14 M o keV/ 

(2.6 + 0.8 (+0.3)) x 10 14 M G , 

and hence M^ Q - / 50 o^200 ■ M 500 - (3.6 ±1.1 (+0.4)) x 10 14 M o 
(ME 5). 

3.2. Optical properties 

We now turn to the optical cluster properties as derived from the 
VLT/FORS2 imaging and spectroscopic data and the wide-field 
multi-band LBT imaging observations. 

The rich galaxy population of XMMU J 1230.3 +1339 is easy 
to recognize when inspecting the top panels of Fig.Q] either as 
a galaxy overdensity measure in the single band (left panel), or 
in color space (right), where the red early-type cluster galaxies 
clearly stand out from the foreground/background. An interest- 
ing note is the fact that the angular size of the foreground Virgo 
galaxy NGC4477 is almost identical to the cluster total radius. 
At a projected cluster-centric distance of about 700 kpc towards 
the West, the optical light of the NGC4477 halo starts to sig- 
nificantly contaminate the colors of background galaxies, and at 
>850 kpc in this Western direction the region is basically masked 
out by the foreground elliptical. The proximity of the cluster 
to NGC4477 hence has two effects on the optical analysis: (i) 
galaxy counts have to be corrected for the geometrically covered 
area in the outer regions, and (ii) the measured colors of galaxies 
overlapping with the outer halo of NGC 4477 are likely biased 
towards the blue. 

Throughout this Sect. 13. 21 we are facing the task to appro- 
priately re-scale measurement apertures that were defined in the 
local Universe to be applicable at z ~ 1 . This is not a straight- 
forward exercise, since (i) one has to distinguish fixed-size phys- 
ical apertures from relative cluster apertures that evolve with 
redshift, and (ii) this evolution depends on the reference den- 
sity relative to which all cluster quantities are defined. For the 
overdensity definition with respect to the critical density p CI (z), 
used in this work, the radius scales in the sel f-similar mo del for 
a cluster at fixed mass as r oc i.e. it 

changes by a factor of 1 .44 from z = to the cluster redshift. On 
the other hand, a definition relative to the mean density p medn (z) 
implies a scaling of r oc (1 + z )~ l , i.e. a decrease by 1.98, and 
with reference to the top-hat collapse density A(z) radii scale as 



r oc E~ 2 ^{z) ■ A~'/ 3 (z), equivalent to a shrinking by about a fac- 
tor of 1.67 at the cluster redshift. For the following discussion 
we c hoose an intermediate radius re-scaling approach consistent 
with lCarlberg et af] d 19971) rocE^iz), i.e. where appropriate we 
apply a radius correction factor of /|- esc = l/E{z) - 0.58. Since 
the total enclosed mass close to the outer radius scales approx- 
imately as M tot (< r) oc r (see e.g. Equ.0, the use of different 
re-scaling approaches for the measurement aperture of optical 
quantities can lead to systematic total mass shifts of +15-20%. 
Working with fixed physical apertures, on the other hand, could 
result in an overestimation of the total cluster mass by up to a 
factor of 2 at z ~ 1, when using optical observables such as rich- 
ness or total luminosity. 



3.2.1. Richness measures 

As starting point, we determine t he classical optical richness 
measure, the Abell richness R 1 ich (lAbelllll958h . defined as the 
number of cluster galaxies within the projected Abell radius 
Ra — 2. 14 /z~ Mpc and with apparent magnitudes in the in- 
terval [mi, mi, + 2], where is the third brightest cluster 
galaxy. We use the LBT/LBC z'-band as the reddest opti- 
cal wide-field imaging filter, providing the highest cluster-to- 
foreground/background contrast and allowing background esti- 
mations in the same field. The magnitude interval for the Abell 
richness in our case is [m^ -, m^z + 2]=[20.74, 22.74], and for the 
re-scaled Abell radius we assume /?A,resc — fresc ' #a - 1240 kpc, 
corresponding to an angular scale of 156". 

We determined the statistical background/foreground galaxy 
counts in this magnitude interval in 17 external background re- 
gions with a radius equivalent to /?5oo each, yielding a median 
z'-band field galaxy density of 8.4 ±1.5 arcmin~ 2 . For the clus- 
ter region within R500, which is not influenced by NGC 4477, 
we obtain a galaxy overdensity of (74 ± 9). The richness esti- 
mate for the scaled Abell radius #A,resc is then (109 ± 17), after 
applying a 10% correction (~3 galaxies) for the geometrically 
covered area in the outer ring beyond /?5oo- This implies that 
XMMUJ1230.3+1339 falls in the Abell richness class R lich = 2, 
defined as having 80-129 counted galaxy members, and would 
thus rank among the top 20% of the richest clusters even in the 
local Abell catalog. 

A more modern and b e tter calibrated rich ness measure is 
N 2 oo dKoester et alJ l2007bt iReves et all l2008b . the number of 
E/S0 ridgeline member galaxies within /?200> fainter than the 
BCG and brighter than 0.4 L*. Thi s richness definition h as been 
used for the MaxBCG catalog of iKoester et alJ (l2007al) of op- 
tically selected groups and clusters from the Sloan Digital Sky 
Survey (SPSS) in the redshift range 0.1 < z < 0.3. We follow 
IReves et aO d2008l) for a self-consistent procedure to determine 
N200, which has been used and calibrated with a different defini- 
tion of #200' eva l uate d relative to the mean density p m ean of the 
Universe. For the selection of E/S0 ridgeline galaxies we use the 
color cuts R-z=2.05 + 0.2 as shown in Fig. [2 with a compara- 
ble color width as used for the MaxBCG catalog. The consid- 
ered z-band magnitude range is then [19.9 mag, 22.9 mag], con- 
fined by the measured BCG magnitude and a SSP (Simple Stellar 
Population) model magnitude of a 0.4 L* galaxy at the cluster 
redshift (see below). For the iterative process, we first determine 
the number of galaxies fulfilling the color and magnitude cuts 
within a fixed projected radius of 1 hr x Mpc ^ 1 .43 Mpc yielding 
N\ h-i Mpc = 56. From the very similar scaled final search radius 
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R* 2Q0 - 0.223 A^-Zrescfe^Mpc - 1.44 MpcQ we obtain the 
final richness measure A^oo - 57, which again includes a small 
correction (3 galaxies) for the region covered by NGC 4477. 

A total mass estimate based on this richness estimate is 
then obtained from the calibrated MaxBCG scaling relation of 
iReves et all (120081) . Due to different radius definitions, the clus- 
ter mass M' 2QQ resulting from this scaling relation is systemati- 
cally higher than M200 used in this work. In the local Universe 
M' 20Q is about 40% larger, which reduces to approximately 8% 
at z ~ 1, due to the fact that p mea n(z) converges to p a -(z) with in- 
creasing redshift. For consistency with our M200 mass definition 
we apply a correction factor /™ ean (z ~ 1) — 0.92 and convert the 
mass units to Ai'Mq 

,aj ,1.16±0.09 

«a» = /cT n ■ (2-03 ± 0.1 1) • (-g£) x 1O 14 M (7) 

= (6.3+ 0.7) x 10 14 M© . 

In terms of systematic uncertainties, we consider an additional 
±15% error for the aperture re-scaling, and +20% for the se- 
lected color cuts and possible foreground/background contami- 
nation in th e richness estimate . Furthermore, all optical scaling 
relations in IReves et aD (120081) are derived from stacked cluster 
quantities, i.e. they do not contain the information on the scat- 
ter of individual clusters with respect to the best fitting relation. 
To account for this intrinsic scatter, we add a conservative min- 
imum of +20% systematic mass unc ertainty to all optical mass 
proxies (see e.g. lPopesso et alj|2005l) . With a combined possible 
systematic mass offset of +32%, the final richness based total 
mass estimate is hence - (6.3 + 0.7 (+2.0)) x 10 14 M o 

(ME 6). 

3.2.2. Brightest Cluster Galaxy 

We address the characteristics of the Brightest Cluster Galaxy 
(BCG) of XMMU J1230.3+1339 and compare them to predic- 
tions from Simple Stellar Population (SSP) models. The BCG 
can be readily identified in the color-magnitude-diagram in 
Fig. [2] as the galaxy with total z-band magnitude zf GG = 19.90 
and a color of R-z^2.04. In the lower left panel of Fig. Q] the 
BCG is marked by a green dashed circle at a projected off- 
set from the nominal X-ray center of about 18" correspond- 
ing to 140 kpc. Furthermore, the BCG (ID02 in Tabled seems 
not to be at rest in the cluster reference system, but rather ex- 
hibits a velocity of about +600km/s with respect to the me- 
dian redshift. The lower panels of Fig.Q] also reveal that the 
BCG has a very close companion galaxy with a consistent color 
in the projected view, which could be an indication for on- 
going or upcoming merging activity. We can summarize, that 
XMMU J1230.3+1339 has a dominant BCG in the cluster core 
with properties in the CMD (magnitude, color, luminosity gap) 
very similar to most local clusters, but additional characteristics 
(cluster center offset, cluster velocity offset, close companion) 
that point towards a dynamical BCG state which is yet to reach 
a final equilibrium configuration at the minimum of the cluster 
potential well. 

Interestingly, we can identify a second galaxy with BCG 
properties on the outskirts of the cluster, from now on referred to 
as BCG2 for distinction. This galaxy can be seen in the Southern 

1 1 Note that the optically defined radius R^g for the determination of 
N200 an d L 2 oo is different from the T x derived radius R 2 oo used through- 
out most of the paper. 



half of Fig. Q] (top right panel) and in Fig.QT] The BCG2 spec- 
trum is shown in the lower panel of Fig.|4] and further proper- 
ties are specified in Table|2] (ID08). With a total magnitude of 
z^ G2 = 20.1, the BCG2 is only 0.2 mag fainter than the core 
BCG and has a slightly bluer color in comparison. Due to the 
projected distance of about 70" - 550 kpc from the cluster center, 
the BCG2 appears in the CMD of Fig. [2] as a spectroscopically 
confirmed (open circle) non-core member (small black dot). The 
cluster restframe velocity of the BCG2 is consistent with zero, 
giving rise to the conclusion that the BCG2 has just entered the 
Rsoq region of the cluster approximately along the plane of the 
sky, which will be further discussed in Sect. 14.2.21 

The comparison to stellar population synthesis models al- 
lows us to probe in more detail the BCG properties and 
to compare them with observations in the local Universe. 
We have computed a grid of S SP models using PEGASE2 
dFioc & Rocca-Volmerangell 19971) with different formation red- 
shifts and metallicities for relating them to the observed red- 
shift evolution o f passively evolving galaxies (for details, see 
lFassbenderil2007l) . The observed R-z color is well matched by a 
model with a single star formation burst at zj — 5 and slightly su- 
persolar metallicity for the BCG and about solar metallicity for 
BCG2. For z~ 1, our models further predict an apparent magni- 
tude of z — 2 1.9 mag for a passive L* galaxy, implying that the 
BCG has a luminosity of 6.5 L* (or m*-2). 

For the comparison of passive galaxy luminosities at higher- 
z to local measurements such as MaxBCG, two different ap- 
proaches could in principle be followed. The first option is to 
perform the measurements in the same observed band, r in the 
case of MaxBCG, and then account for passive evolution and a 
K-correction term of order 2 mag using the model predictions. 
The observationally preferable second approach uses the best 
suited band for the measurement, z in our case, then applies sig- 
nificantly smaller evolution and K-correction terms, and finally 
transforms the local values into the desired reference band using 
the model. Here we follow the latter method in order to take ad- 
vantage of the ~2 magnitudes brighter passive cluster galaxies 
in the z-band compared to r or R. 

We determine the absolute z-magnitude by applying 
the distance modulus equation tobcg(z) - M^cg = 25 + 
5 log (d lma [Mpc]) + K BCG (z)- With zf cc = 19.90, the luminosity 
distance d\ um = 6403.6 Mpc, and a K-correction of Kbcg — 0.82 
we obtain M^ CQ = -24.95 + 0.03 mag for our assumed con- 
cordance co smology. Using the absolute magnitude of the sun 
Mq = 4.51 dBlanton et alj|2003l) we can derive the BCG lumi- 
nosity as i BCG - 6.1 x 10 11 Lq. 

The BCG luminosity is know n to correlate weakly with the 
cluster mass dLin & Mohrll2004j). Althoug h found to be a fairly 
poor tracer of M200, Reyes et ail d2008l) provide a calibrated 
Lbcg-M scaling relation for BCG r-band luminosities evaluated 
and K-corrected at z = 0.25. With the discussed correction term 
ymean an( j i um j nos ities and masses transformed to Lq and 

Ii^q Mq we arrive at the following relation 

/ jjr n 1.10±0.13 

<o CG = / c T an (l-53±O.lO)[ 1O2x B 1 C g llL0 j xlO 14 M .(8) 

In order to apply this relation, we need to passively evolve the 
BCG to z = 0.25 and then evaluate its r-band luminosity (Mq = 

4.76) yielding U BCC =* 1.7 x 10 11 L©. From Equ.[8]we obtain the 
mass estimate M 20 o -(2.4 + 0.4) x 10 14 M© (ME 7a). 

So far, we assumed that the BCG's stellar population is 
passively dimming without a change in stellar mass. If we in- 



R. Fassbenderetal.: A pan-chromatic view of XMMU J1230.3+1339 at z = 0.975 



13 



stead assume a non-evolving absolute magnitude, i.e. the dim- 
ming of stellar light is compensated by stellar mass acquisi- 
tion, then the r-band luminosity at z — 0.25 is 70% higher 
(Lg CG 2.9 x 10 11 L0) resulting in a mass estimate of M200 — 
4.4 ± 0.7) x 10 14 M© (ME 7b). For consistency with lReves et al.l 
20081 we use the non-evolving luminosity as reference mass 
proxy, and treat the offset to the result including passive evo- 
lution as a systematic error. With an additional +10% system- 
atic for the SSP model evolution uncertainty and +20% for the 
intrinsic scatter, we arrive at the final cluster mass estimate of 
M 200 G ^( 4 - 4 ± O.7(!^))xlO 14 M (ME 7). 

3.2.3. Optical luminosity 

The next optical observable and useful mass proxy is the clus- 
ter luminosity L20O1 defined as the summed r-band luminosity 
of all red ridgelin e galaxies in ./V200, K-corrected to z = 0.25 
dReves et al.l 12008). Following the procedures of the last two 
Sects. [22jJ and 13.2.21 we derive a total ridgeline z-band lu- 
minosity inside R* 2QQ of L z m =* 9.2 x 10 12 Lq corresponding to 
15.1 Lg CG . In contrast to the BCG, where the accretion of addi- 
tional stellar mass can potentially compensate the dimming of its 
predominantly passive stellar population, the natural choice for 
the treatment of the total ridgeline light, which should be closely 
linked to the total stellar mass in the cluster, is the assumption of 
passively evolving luminosities. The expected evolved and trans- 
formed bulk luminosity in the r-band at z = 0.25 is therefore 
L!j 00 ^2.5xl0 12 Lq. The corresponding scaling relation (in units 
of L© and M©) 

/ tt \1.40±0.19 

=fc7™(2-5 1*0.24) | g 16x 2 1 °° llL0 ) xl0 14 M©(9) 

yields a total cluster mass estimate of M 2 oo -(11.1 ±2.8)x10 I4 Mq 
(ME 8a). The systematic uncertainties now add up to about 
±34%, based on the discussed contributions from aperture re- 
scaling (15%), color cut selection (20%), the SSP evolution 
model (10%), and the intrinsic scatter assumption (20%). The 
final L200 mass proxy for the cluster is hence (11.1 ± 

2.8(+3.7))xl0 14 M© (ME 8). 

While the passively evolved BCG luminosity pointed to- 
wards a low cluster mass estimate (ME 7a), the total L200 light 
(ME 8) would give rise to the opposite conclusion of an ex- 
tremely high total mass. This discrepancy is to be attributed to 
the highly populated bright end of the red galaxy ridgeline (see 
Fig.0, which apparently places XM MU J1 230.3+1339 in a low 
mass-to-light ratio regime (see Sect. l3.2.5T > with respect to the 
underlying MaxBCG cluster sample used for the calibration of 
the L200 relation. 

3.2.4. Optimal optical mass tracer 

iReves etail d2008l) also tested the possibility of an improved op- 
tical mass tracer by considering power law combinations of ei- 
ther the richness N200 or the total luminosity L200 with the BCG 
luminosity Lbcg- The resulting optimal optical mass tracer with 
an improved scatter of the form (in units of tc} L© and hzi M©) 

,aj >1.20±0.09 

M2°oo=/c7 an (l-81±0.11)(^) (10) 

/ , v 0.71 ±0.14 

■ _ Lbcg x10 14 M© 

\LbCg(N20o)/ 



points towards an anti-correlation of the richness with the 
BCG luminosity at fixed cluster mass. Here Z,bcg(^V20o) is the 
mean local BCG luminosity at a given richness Lbcg(A^20o) = 
3.14 • Moo' x 10 10 L©. Using the determined parameters for 
XMMUJ1230.3+1339 with the non-evolved BCG luminosity 
Lg CG 2.9 x 10 11 L© and N200 - 57 we obtain an optimal op- 
tical mass estimate of M^(S.6 + 1.3 (^))xl0 14 M© (ME 9). 
The systematic uncertainties here include the richness system- 
atics of Sect. l3.2.T1 and the passive evolution offset discussed in 
Sect. l3.2.2l which would decrease the total mass by 32%. 

3.2.5. Mass-to-light ratio 

With our available LBT wide-field imaging data, we have an 
independent way to confirm the high optical luminosity of the 
cluster using the statistical background subtraction method. This 
approach does not include a color cut, i.e. also galaxies out- 
side the red ridgeline region are counted in a statistical sense. 
We focus here on the uncontaminated /?soo region of the clus- 
ter and determine the background luminosity in the 17 external 
regions as in Sect. 13. 2. Tl taking all galaxies with m > hibcg into 
account. This way we measure a total background subtracted ap- 
parent z'-band magnitude of = 16.79 mag within the ^500 

aperture, corresponding to a restframe luminosity of L l °^, dp = 

(1.06 + 0.16) x 10 13 L© (see Sect. l3"T2l . The integrated appar- 
ent magnitude of the cluster is comparable to the total fore- 
ground/background light contribution, implying that about half 
of all detected z'-band photons in the R500 region originate from 
XMMU J1230.3+1339. Only about 6% of this luminosity is con- 

tributed by faint galaxies with magnitudes m z > 22.9. The red 
ridgeline galaxies of Sect. 13. 2. 31 account for approximately 61% 

r 

of the galaxy light with m z < 22.9 in the R500 region. 

Although the system is at z ~ 1, the contrast of the clus- 
ter's z'-band light with respect to the integrated background is 
sufficiently high to allow an evaluation of the spatial light dis- 
tribution in the core region and some surrounding substructures 
without any particular pre-selection of member galaxies. As be- 
fore, we merely mask out stars and galaxies brighter than the 
BCG (ot<otbcg) in the sky-subtracted z'-band image and apply 
Gaussian smoothing with a 180kpc (22") kernel over the field. 
The resulting observed z'-band light distribution of the cluster 
field is shown in the lower left panel of Fig. [9] displayed with 
6 linear spaced contours corresponding to significance levels of 
2cr-12cr above the background. 

Next, we can derive the total z'-band mass-to-light ratio 
within the R500 aperture using the final best mass estimate of 
the cluster (see Sect. l5.2b 

^) » = (46.7 ± 11.3)^. (11) 

Lz >™w / 5 oo^2oo • /.:;;;-,„ L ° 

Here /3d^ a p = A^soo.ap/^oo - 1-65 is the ratio of the aperture 
mass M500 ,ap within a cylinder with radius 7?goo and the three 
dimensional spherically en closed ma s s M500 , following the pro- 
jection formula provided in lJee et a! 1 (120051) for the derived sin- 
gle y6-model parameters of the cluster. In order to compare this 
rather small value with locally determined mass-to-light ratios 
of clusters, we have to take evolution effects into account, since 
the optical luminosity is expected to undergo significant changes 
over a time span of 7.6 Gyr due to the increasing mean age of 
the underlying stellar populations. A suitable approximation for 
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this change in the integrated luminosity is the passive evolution 
model, since the total light is dominated by galaxies on the red 
ridgeline or very close to it in color space. This model predicts a 
dimming of the absolute z'-magnitude by 0.7 mag, equivalent to 
a decrease of the luminosity by a factor of 1 .9 at z — 0, resulting 

in an evolved mass-to-light-ratio of (M/L z ap = 89 (Mq/Lq). 
The transformation into the V-band leads to an additional in- 
crease to (^/i V )5oo ap - 184(M©/L©), i.e. to a mass-to-light 
ratio in the lower half of the Mj L p arameter region span ned by 
other massive local clusters (see e.g. lPopesso et alj|2007l) . 



3.2.6. Velocity dispersion 

Under the assumption of virial equilibrium a dynamical mass 
estimate can be obtained from a measurement of the cluster's 
radial velocity dispersion cr,. (e.g. iFinn et all 120051: iKurk et all 
120091) 



M 2 oo 



3 a 2 R200 1 .7 



O r 



— (- 

E(z) UOOOkms 



T ) 3 xl0 15 /7 7 - 1 M o . (12) 



3.2.7. Spatial galaxy distribution 

As the final optical result discussed in this work we derive a spa- 
tial density map of the red ridgeline galaxies selected with the 
previously applied R-z color cut as shown in Fig. [2] The density 
map was constructed by marking the world coordinate positions 
of the 88 selected galaxies over the FORS 2 field-of-view and 
the subsequent application of an adaptive smoothing filter. The 
resulting spatial red ridgeline galaxy distribution is displayed in 
the upper left panel of Fig.|9]with log spaced color cuts and con- 
tour levels analogous to the representation of the X-ray surface 
brightness in the upper right panel. Going from the lowest to the 
highest contour level corresponds to a galaxy density increase 
by a factor of 26, with the peak density coinciding with the X- 
ray center marked by the white cross. Within the R^m region 
of the cluster, two additional galaxy concentrations towards the 
North- West and South-East from the center are obvious. Along 
the same axis, galaxy overdensities can be traced out to beyond 
the nominal cluster radius. These components and features will 
be further discussed in the multi-wavelength view of the cluster 
in Sect.g] 



We determined the velocity dispersion for 
XMMUJ1230.3+1339 us ing the 'ro bust' estimator of 
iBeers et al.l d 19901) and iGirardi et alj d 1993b applied to 
the 13 secure cluster members! 12 l in Table|2] and found 
o> ^ 658 + 277 km s -1 (Fig.|5]> in good ag r eemen t with 
the classical dispersion measure of iDanese et al.l (fl980T) . From 
Equ.[T2] we derive a formal dynamical cluster mass estimate of 
M 200 ^ 2 - 8 -23 x 10 14 M Q (ME 10). 

Although the statistics of available galaxy member redshifts 
is low with corresponding large uncertainties, the general sit- 
uation for a first dynamical mass estimate based on about a 
dozen redshifts is typical for most optical/infrared selected dis- 
tant galaxy clusters (e.g. IStanford et al. 120051: lEisenhardt et alJ 



I2008t IWilson et al.l 120091: iMuzzin et a" 



2009). While at local 



redshifts dynamical mass estimators can provide robust mea- 
surements for relax ed clusters even for rather small spectro- 
scopic samples (e.g. iBiviano et al.ll2006l) . high redshift systems 
are particularly prone to be affected by potentially significant bi- 
ases. In the case of XMMU J1230.3+ 1339 the fairly low velocity 
dispersion estimate compared to similarly rich and massive local 
clusters becomes plausible in the light of the large-scale struc- 
ture and component analysis in Sect. 14. 21 More than half of the 
available redshifts are likely to be attributed to galaxies resid- 
ing in coherently infalling structures, i.e. they are not virialized 
tracer particles of the underlying potential well. What is effec- 
tively measured is hence a combination of the 'fingers-of-God 
effect', corresponding to the virialized tracer particles the mass 
estimate is based upon, and the prolate ellipsoid in redshift space 
due t o the coherent infa ll of the surrounding large-scale structure 
(e.g. iGuzzo et al.ll2008h . The effect of this latter systematic ve- 
locity dispersion bias depends on the angle of the targeted bulk 
matter flows with respect to the plane of the sky, which can sup- 
press the measured velocity dispersion in the case of coherently 
infalling substructures mainly in the transversal direction (see 
Fig.[T2l. or boost up the dynamical mass estimates for radial co- 
herent flows. 



3.3. Weak lensing results 

The projected total mass surface density is observationally 
accessible through the method of weak gravitational lens- 
ing. We have used the deep multi-band LBT imaging data 
of Sect. l2.1.2l to perform a full weak lensing analysis in the 
field of XMMU J1230.3+1339, which is presented in detail in 
the accompanying Paper II. For obtaining a complete multi- 
wavelength view on the cluster, we summarize here some rel- 
evant results of the analysis for this work. 

With a redshift of z ~ 1, XMMU J1230.3+1339 is currently 
at the feasibility limit for weak lensing studies using ground- 
based imaging data. However, the available LBT data allowed 
the detection of a significant weak lensing signal based on the 
shape distortions of background galaxies, which were separated 
from foreground objects by means of multi-band photometric 
redshifts. The resulting weak lensing signal-to-noise map, which 
is proportional to the projected total mass surface density, is dis- 
played in the lower right panel of Fig. [9] The five linearly spaced 
contour levels span the significance range from half the peak 
value (1.75 cr) to the maximum at 3.5 cr. The location of this 
mass density maximum coincides with the X-ray center and the 
peak of the galaxy overdensity. Furthermore, the general cluster 
elongation in the SE-NW direction is also clearly visible in the 
weak lensing map. The detected extension towards the North - 
West is significant to better than 2 <x, whereas the South- Western 
feature is to be considered tentatively due to the close proximity 
to the contaminating foreground galaxy NGC 4477. 

For the reconstruction of the total spherical 3D-mass of the 
cluster we fitted two parametric mass models to the observed 
tangential shear profile (see Paper II for details). For an assumed 
singular isothermal sphere (SIS) density profile we determine the 
corresponding best fitting velocity dispersion of the cluster to be 
of s ^ 1271 + 255 km s" 1 . Using R 2 oo of Sect. mTTl as outer clus- 
ter boundary we obtain a SIS modeFl mass estimate (ME 11) 
of 



12 This method yields a slightly lower estimated cluster redshift of 
Zd = 0.9737 compared to the median redshift z = 0.9745 of Sect. l2.1.3l 
but still well within the 1 cr uncertainty. The first value is likely biased 
low due to coherently infalling sub-structure (see Sect. l4.2.ll . so that 
the median value is preferable as cluster-centric reference point. 
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Fig. 9. A pan-chromatic view of XMMU J 1230.3+ 1339 seen through projected surface density maps of different physical quantities 
within the cluster's R200 region (outer dashed circle). All panels show the same 4.8'x4.8' sky region in standard orientation (North is 
up, East to left), the nominal cluster center (central cross) and R500 (inner dashed circle) are also indicated. Top left (a): Density map 
of color-selected red galaxies (log-spaced contours). Top right (b): Chandra X-ray surface brightness map (log-spaced contours). 
Bottom left (c): Total z'-band light distribution map (linearly-spaced contours). Bottom right (d): Weak lensing mass density map 
(linearly-spaced contours). 



* (7.6 ±3.3 (+3.0)) x 10 14 M o . 

As a second parametric model we considered a Na varro-Frenk- 
White (NFW) density profile (iNavarro et alJll997h with a con- 
centration parameter of c = 4.0 and a scale radius of r s = 
354kpcQ These NFW parameters were determined from the 
measured shear components, in good agreement with the ex- 
pected values from simulatio ns of halos of sim i lar ma ss and red- 
shift bv lBullocketai] ( |2001l) and lDolag etall d2004l) . This best 
fit NFW model to the weak lensing data yields an enclosed to- 
tal mass within R 2m of (6.9 + 3.3 (+2.2)) x 10 14 M 
(ME 12). 



For consistency and robustness we used the X-ray deter- 
mined value for R-200 as outer cluster boundary. However, this 
way the weak lensing mass estimate is not fully independent. 
Indeed, the self-consistent direct determination of #200 from the 
best fitting NFW profile yields a significantly larger radius of 
/?2 ( x) W - 1415 kpc. To account for this potential positive offset in 
the outer cluster radius, we considered the resulting mass differ- 
ences as additional systematic uncertainty. 



For hyo, slightly scaled from the values given in Paper II for ^72. 
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Fig. 10. Cluster components and large-scale structure environment of XMMU J1230. 3 + 1339 in a 7.0' x 7.3' LBT color composite 
view with different density contours overlaid. Six identified cluster components are marked by dashed green boxes labelled with 
roman numerals. Left: Contours show the galaxy densities in red and the weak lensing total mass density in white. Right: Same field 
with the z'-band light distribution in yellow and XMM-Newton X-ray contours in cyan. 



4. Multi-wavelength view on large-scale structure 
and small-scale physics 

We now combine the results derived for the different cluster 
components (i.e. galaxies, ICM, and Dark Matter) in order to 
obtain a detailed pan-chromatic view of XMMU J1230.3+1339. 

4.1. Comparison of projected density maps 

Figure|9]displays the four derived spatial density maps (Sect. [3} 
within the cluster's /?2oo region (outer dashed circle), corre- 
sponding to an angular image side length of 4.8'. The maps 
are based on independent data sets or techniques and visualize 
different physical projected properties and matter components. 
The top panels show the log-spaced density contours of the red 
ridgeline galaxies (left, a), based on the VLT/FORS 2 data of 
Sect. l2.1.Tl and the Chandra observed (Sect. 12. 2^21 X-ray sur- 
face brightness distribution (right, b). The bottom panels display 
the linearly-spaced contours of the total projected light distribu- 
tion in the z'-band (Sect. l3.23l l on the left (c) and the weak lens- 
ing signal (Sect. l3.3l ) on the right (d) both derived from the deep 
LBT/LBC imaging data of Sect. l2.1.2l (see Paper II for details). 
In terms of physical projected cluster quantities, the WL signal 
is proportional to the total mass surface density (d), the observed 
z'-band flux represents the restframe B/V total stellar light (c), 
the X-ray surface brightness is a measure of the projected gas 
density squared n 2 e (b), and the galaxies can be approximately 
considered as tracer particles of the underlying gravitational po- 
tential well (a). 

The first observation is that the nominal cluster center (cen- 
tral cross), measured as the 'center-of-mass' of the X-ray emis- 
sion detected with XMM-Newton (Sect. 13. Il l, shows a very good 
correspondence to the peak signal within a few arcseconds in all 



projected density maps. This implies that the galaxies, the ICM 
gas, and the dark matter component are linked to the same (pro- 
jected) center of their corresponding density profiles. Note that 
this concordant cluster center for XMMU J1230.3+1339 does 
not coincide with the BCG location (see lower left panel of 
Fig-tO-. which is often used as fiducial reference center in optical 
cluster studies. 

A second common feature in all maps is an elongation or 
extension in the SE-NW direction, which indicates the appar- 
ent main axis for the cluster assembly and the connection to the 
large-scale structure filaments of the surrounding cosmic web. 
Just NW of the center, all projected density maps (at lower res- 
olution for the WL signal) show another clear extension feature 
in the core region, which is identified with an ongoing merging 
event discussed in Sects. l472l S d4.4l The galaxy (a) and total mat- 
ter overdensities (d) in the NW direction can even be traced to 
beyond the nominal R200 cluster radius, indicative for the pres- 
ence of a LSS filament. In the opposite SE direction, the galaxy 
map (a) reveals a further density peak just inside the inner R500 
radius, also reflected in the light map (c), and another extension 
to beyond the cluster radius. 

Towards the South- Western (SW) direction from the system 
center, an additional but less pronounced likely cluster feature is 
present in all maps. In contrast to the conclusions on the over- 
all shape of XMMU J1230.3+ 1339 along the identified main 
NW-SE axis, which was not not influenced by NGC 4477 (see 
Fig.fTOb. the proximity to the bright foreground galaxy in the 
SW region now hinders the component analysis of the cluster in 
the far background and weakens the significance due to the de- 
teriorated achievable shape and color measurement accuracy of 
background galaxies and possible residual contributions to the 
total light and X-ray emission in this region. As a result, the 
measured South- Western red galaxy density is to be regarded as 
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lower limit (a), while the X-ray emission (b, see also the deeper 
XMM SNR map in Fig.[l0] right) and total light distribution (c) 
could be slightly increased, and a tentative character applies to 
the SW feature in the weak lensing map (d). 

4.2. Cluster components and large-scale structure 

We turn the focus on the global view and the characterization 
of the dynamical state of XMMU J1230.3+1339 based on the 
available multi-wavelength data set. Figure[10] displays an opti- 
cal color composite with a larger 7.0' x7.3' field-of-view and 
the different projected density contours overlaid. The weak lens- 
ing signal (white) and the galaxy densities (red) are shown in 
the left panel, and the z'-band light distribution (yellow) and the 
XMM-Newton X-ray emission (cyan) on the right. For a bet- 
ter qualitative evaluation of low surface brightness X-ray emis- 
sion in the cluster outskirts the combined XMM-Newton data 
is used (Sect. l2.2TTb resulting in an improved sensitivity at the 
cost of lower spatial resolution compared to Chandra. The low- 
est four cyan contours correspond to significance levels of the X- 
ray emission of (0.5, 1, 2, 3)cr above the background, whereas 
the central cluster contour displays a 35 cr significance. 

We have identified a total of six different cluster components 
or associated structures in the multi-wavelength maps, with iden- 
tification labels I- VI in order of decreasing significance: the clus- 
ter core (I), an ongoing merger event towards the NW of the cen- 
ter (II), two infalling groups on the cluster outskirts (III & VI), 
and two filament extensions beyond the nominal cluster radius 
towards the South-East (IV) and the North- West (V). 

The central cluster region, indicated by box I in Fig.[10j en- 
compasses the projected density peaks of all the main matter 
components, i.e. the Dark Matter dominated total mass density 
(white contours), the hot ICM gas component (cyan), and the 
cold baryon component in form of the total stellar light (yellow) 
or the spatial galaxy density (red). The region also includes the 
spectroscopically confirmed BCG (ID 02 in Tabled at a radial 
distance of 18" from the nominal cluster center (bottom panels 
of Fig.[T]i. For the further discussion, we follow the flow of ac- 
creted matter onto the main cluster halo from the largest scales 
towards the central regions. 

4.2.1 . Cosmic web filaments 

The North- Western galaxy filament (box V) is traced to a ra- 
dial distance of more than twice R200, where the current limit 
is set by the FORS 2 FoV, i.e. the actual overdensity of red 
galaxies might extend even further. The projected galaxy den- 
sity appears to have a minimum close to R-200 and then increases 
again towards the local maximum in the filament of region (V) 
(left panel of Fig.lT0T>. Three spectroscopic cluster members (IDs 
09,10,12 in Tabled are located in the cluster-filament transition 
zone, all are blue shifted with respect to the median redshift with 
cluster-centric velocities in the range [-800,-1 300]km/s, which 
is indicative for a line-of-sight bulk flow velocity of the filament 
structure of -1 OOOkm/s. In terms of the weak lensing mass 
density, the signal of an extended mass distribution is measured 
well into the cluster-filament transition regime at a signal-to- 
noise level of more than two. From the X-ray side, there are some 
indications for low surface brightness extensions from the clus- 
ter core towards the NW filament (right panel of Fig.lTOb. and 
several X-ray point sources in box V, which may or may not be 
related to the cluster environment. 




Fig. 11. Closeup view on the infalling BCG2 group (III). 

The opposite galaxy extension towards the South-East (IV) 
is less pronounced and reaches about 0.5 Mpc beyond /?200- The 
association of the three X-ray sources with corresponding light 
overdensities in box IV to the filament can currently not be 
probed with the limited spectroscopy available, which includes 
one secure member galaxy (ID 1 1 in Tabled matching the me- 
dian redshift, i.e. with vanishing cluster-centric line-of-sight ve- 
locity. The tentative interpretation of a galaxy filament along the 
plane of the sky is fostered by the similar cluster-centric veloc- 
ity of BCG2 (ID 08) as part of the infalling galaxy group III just 
inside R500. 

4.2.2. Infalling groups on the cluster outskirts 

BCG2 (see Sect. 13.2.21 appears to be the central galaxy of an in- 
tact group environment approaching the cluster center from the 
SE. FigureQTj displays a closeup view on this group component 
(III), which still seems to be a bound sub-system without signs 
of strong interactions with the main cluster components yet. The 
very low cluster-centric velocity of BCG2 points at a radial in- 
fall orbit towards the center of XMMU J1230.3+1339 along the 
plane of the sky with a current cluster-centric distance of about 
560 kpc. The BCG2 group is clearly visible as a local peak in 
the red galaxy density (a) or total light (c) map of Fig. [9] and a 
less pronounced distortion in the XMM-Newton X-ray contours 
(right panel of Fig.fTOb. At the current mass and spatial resolu- 
tion of the weak lensing map, the BCG2 group does not produce 
a local mass peak feature in the weak lensing contours, suggest- 
ing that the associated total mass of the BCG2 group is small 
compared to the main cluster halo. 

A second infalling group candidate at a slightly larger pro- 
jected cluster-centric distance of about 900 kpc is marked as 
component VI in Fig. [10] This South- Western cluster component 
is currently the least secure one, as spectroscopic confirmation 
is still lacking and the discussed proximity to NGC 4477 may 
impose a bias on the projected quantities. Nevertheless, all four 
density maps show features that point towards an off-axis infall 
scenario of matter also from the SW direction, but currently un- 
constrained along the radial direction. Besides the tentative mass 
feature in the weak lensing map (white contours) in the direction 
of component VI, XMM-Newton detects X-ray emission within 
this region on the 3 cr significance level (cyan). This X-ray emis- 
sion seems to be extended, but due to the low-surface brightness 
a further characterization is currently not feasible. The general, 
almost tangential direction with respect to the cluster center is 
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also outlined by the coincident total light distribution of this 
region (yellow) and the density of red galaxies (red), which is 
certainly biased low compared to the true density as a result of 
color contamination from the blue foreground emission on the 
outskirts ofNGC4477. 

4.2.3. Central merging activity 

Moving back into the core region of XMMU J1230.3+1339, we 
now turn to the cluster component II, which extends to the North- 
West of the nominal center at a projected cluster-centric dis- 
tance of 100-400 kpc. This component features prominently in 
the galaxy density (a) and total light (c) maps of Fig.[9]as linear 
NW prolongation close to the central peak value. The Chandra 
X-ray surface brightness map (b) reveals an associated wedge- 
like region with an almost constant SB, followed by a rapid out- 
ward SB decline at a projected distance of about 250 kpc. This 
surface brightness plateau was already visible in the azimuthally 
averaged radial SB profile of Fig. [8] and can also be recognized 
at lower spatial resolution in the XMM-Wwtofi X-ray contours 
of Fig. [10] The weak lensing map (d) provides a view on the 
total mass distribution at a further reduction of resolution, but 
still reveals a significant elongation of the projected mass in the 
direction of component II. 

The overall geometry and features in the projected den- 
sity maps are similar to the merger configuration in the well 
studied 'Bullet Cluster' IE 0657-56 at a lower redshift of z — 
0.3 (e .g. lMarkevitch et al.ll2002tlClowe et all2006tlBradac et al l 
In comparison, the observable signs of merging activ- 
ity in XMMU J1230. 3+1339 occur at smaller angular scales of 
551' and correspondingly shorter projected physical distances of 
$500 kpc, which implies either a steeper viewing angle with re- 
spect to the plane-of-the sky or a merger event closer to cluster 
core passage. 

All currently available redshifts within component II are 
blue-shifted with respect to the median zeropoint. Besides one 
secure member (ID 01 in Tabled at cluster-centric line-of-sight 
velocity of =* -400 km/s, four more tentative redshifts were iden- 
tified closer towards the NW front at larger blue-shifted values 
of about [-700,- 1300,- 1300,-2800] km/s, shown as squares in the 
lower left quadrant of Fig. [5] With the currently available spec- 
troscopy, we can establish a negative line-of-sight bulk velocity 
of the central merger component II, with a tentative best esti- 
mate for the center-of-mass velocity of the merging group of 
v, (-800 + 500) km/s. 

4.3. 3D mass accretion geometry 

Using the identified cluster components in the 2D multi- 
wavelength maps of Sect. 14.21 and the available spectroscopic 
information, we can reconstruct a tentative model of the three- 
dimensional accretion geometry of the cluster environment of 
XMMU J1230.3+ 1339. In Fig.[T2] the horizontal directions cor- 
responds to a cut along the cluster's principal SE-NW axis, and 
the vertical direction indicates the tentatively identified radial 
line-of-sight configuration consistent with the currently avail- 
able limited spectroscopic information for the individual com- 
ponents I-V. Component VI is not considered here since the ori- 
entation along the redshift space direction is not constrained at 
the moment. 

The cluster core (component I) with the nominal center 
(cross) and the offset BCG (black dot) was most likely passed 
by the group bullet (II) within the last lOOMyr at a small but 
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Fig. 12. Sketch of the reconstructed 3D accretion geometry for 
the cluster environment. The horizontal direction shows a cut 
through the main SE-NW cluster axis, whereas the tentative ra- 
dial (i.e. redshift space) distances are plotted along the vertical 
axis for components I-V. The nominal cluster center is indicated 
by the central cross, the positions of the BCG and BCG2 by dots, 
and the reference radii R200 and /?5oo by the solid and dashed cir- 
cles. The physical scale in the cluster reference frame is given 
in the upper left corner and the value in the lower part refers to 
the comoving distance to the cluster in the assumed cosmology. 
Green arrows indicate schematically the approximate traversed 
distance within 200 Myr for linearly extrapolated trajectories un- 
der the assumption of matter infall velocities of 2500 km/s for 
the central merging group II (see Sect. 14.41 ) and fiducial refer- 
ence values of 1000 km/s for the outer components (III,IV,V). 



non-zero impact parameter under a relatively small angle with 
respect to the plane-of-the-sky (see Sect l4.4l ). Judging from the 
undistorted galaxies in the immediate vicinity of the nominal 
cluster center, and the concordant peak location in the density 
maps of the galaxy-, light-, ICM-, and total mass distribution, 
the indications do not support the scenario of a direct hit (i.e. 
zero impact parameter) of the seemingly intact cluster center. 
The green arrow connected to component II visualizes the ap- 
proximate traversed distance of an assumed linearly extrapolated 
bullet trajectory at a relative group velocity of 2500 km/s within 
200 Myr. This distance of ^0.5 Mpc emphasizes the relatively 
short timescale of ~ 100 Myr of a fierce interaction of the clus- 
ter core with the infalling group, and the implied rarity of these 
events to be 'caught-in-the-act' . 

The BCG2 group (component III) has recently entered the 
^500 region of the main cluster from the South-East and does 
not yet show visual signs of interactions with the main halo. The 
spectroscopic information indicates a radial orbit with a mini- 
mal impact parameter. Assuming a reasonable fiducial current 
infall velocity of 1000 km/s, we obtain a traversed distance of 
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about 200 kpc in 200 Myr (green arrow), and a likely direct core 
impact within a timescale of roughly ~0.5 Gyr. For simil ar as- 
sumed fiducial bulk flow velocities (e.g. lPivato et al. 2006) from 
the identified galaxy filaments IV and V, the timescales for mat- 
ter infall to the cluster core are approximately 1 -2 Gyr. From this 
first model of the 3D configuration of XMMU J1230.3+1339 
and its large-scale structure environment, the cluster is likely to 
be in a continuing active mass accretion and growth phase for 
several Giga years to come. 

4.4. The central group 'bullet' 

Merging events constitute a crucial role in the galaxy cluster 
growth history within the hierarchical structure formation frame- 
work. However, the in situ capture of merging events close 
to the cluster core passage is rare due to the relatively short 
core interaction and crossing timescale of the order of 100 Myr. 
Since merger events are important laboratories for a rich num- 
ber of further detailed investigations on Dark Matter proper- 
ties, shock physics of the ICM, galaxy interactions and disrup- 
tions, and relativistic particle acceleration, we summarize the 
current observational status of the new high-z group 'bullet' 
in XMMU J1230.3+1339 here in order to establish a tentative 
multi-wavelength reference model of the merger configuration, 
which is to be tested and extended with upcoming deeper obser- 
vations. 

FigureQj] shows the currently most detailed optical view on 
the dynamically active region of interest to the North- West of 
the nominal cluster center, with Chandra X-ray contours and the 
position of a VLA/FIRST radio source overlaid. 

4.4.1. Merging features in the ICM 

The merging process of a sub-cluster with the main halo is ex- 
pected to occur at peak velocities corresponding to typical Mach 
numbers of M = v/c s - 1-3, where c s is the velocity of sound in 
the main cl uster and v is the relative veloc ity of the colliding gas 
fronts (e.g. IMarkevitch & Vildilininll2007l) . Interestingly, the ex- 
pected typical value is obtained, when interpreting the opening 
angle with respect to the symmetry axis of the observed wedge- 
like structure in the X-ray surface brightness map as Mach cone. 
The dashed yellow lines in Fig.[13]represent the approximate ob- 
served wedges of constant surface brightness with an opening 
angle of <p 28°, which yields M — 1/ sin(i^) 2.1. The cur- 
rently available Chandra data does not have sufficient depth for 
detailed diagnostics of the wedge-like X-ray surface brightness 
discontinuity, which would require measurements of the gas den- 
sity and temperature jumps across the front. However, using this 
typical Mach number of M = 2. 1 as a first estimate for the relative 
velocities, we can obtain a first approximation of the orientation 
angle with respect to the plane of the sky. From the ICM tem- 
perature Tx - 5 .3 keV - 6 .2 x 1 7 K we derive a sound velocity of 
c s =* 1480 • (77 10 8 K) 1/2 km/s = 1160km/s, and hence a relative 
velocity between merging group and cluster of v re i - 2400 km/s. 
The combination with the radial line-of-sight velocity compo- 
nent estimate of v r - -800 km/s yields a first approximation for 
the angle of the merger axis with respect to the plane-of-the-sky 
of ^20°. 

In analogy to the 'Bullet Cluster' IE 0657-56, the observed 
wedge structure in XMMU J1230.3+1339 is likely a 'cold front' 
in the 'stripping stage', i.e. a contact discontinuity on the bound- 
ary of the colder gas of the merging group and the hotter clus- 
ter ICM expected to form close to the core passage point (e.g. 



iMarkevitch & Vikhlininl [20071) . The associated bow shock, as 
seen in IE 0657-56, is expected to be detectable only on the 
outbound trajectory, i.e. it might not have formed prominently 
yet at the particular merger stage in XMMU J1230.3+1339. The 
tentative interpretation of a bullet-like merging event close to 
core passage with an observable 'cold front' (i.e. with an out- 
ward rising temperature) would also be consistent with the in- 
dications of an increasing ICM temperature at radii beyond 
the SB wedge (Sect. l2.2TT1 ). Moreover, in the early 'stripping 
stages' the SB wedge opening angle could be closer connected 
to the actual Mach cone, which is not the case any more for 
cluster IE 0657-56 due to the decreasing bullet size with time 
and a corresponding increase in the apparent opening angle 
( IMarkevitch & Vikhlininl [20071) . A more detailed X-ray analy- 
sis of the merging physics in XMMU J1230. 3+1339 based on 
deeper Chandra observations will be presented in a forthcoming 
paper. 

4.4.2. Interactions of the galaxy component 

Indications for seeing a snapshot of an early merger stage close 
to core passage also arise from a closer look at the galaxy com- 
ponent in Fig. [13] The front end of the X-ray wedge exhibits a 
surprisingly coincident alignment with a dense cloud of galaxies 
(upper right quarter of Fig.[T3l). seemingly the leading front of 
the merging group's galaxy content. The practically collisionless 
galaxy component has thus not yet been spatially separated from 
the fluid-like ICM component, as in the case of IE 0657-56. 

Several galaxies in the central part of Fig. [T3l show signs of 
ongoing galaxy transformation processes. The object with the 
bright core (label i) at the inner tip of the wedge North of the 
yellow foreground spiral (label f) is a candidate event for the 
collision of two gas rich galaxies, based on the strong distor- 
tions of the inner isophotes and the blue color. The photometric 
redshift z p h t - 1.0 of Paper II for this galaxy indicates a tenta- 
tive cluster membership, which is to be spectroscopically con- 
firmed. A different process is seen at work in the red spectro- 
scopic cluster member (ID01 in Tabled to the East of the green 
cross, which features a trail extending about 40 kpc to the SE (la- 
bel t). This galaxy trail could be tidally stripped stellar material 
resulting from the gravitational interaction with the core of the 
mai n halo and could en d up as part of the Intracluster Light (ICL) 
(e.g. lZibetti etalj2005l) . More detailed studies on galaxies trans- 
formation processes in the dynamically active core region of 
XMMU J1230. 3+1339 will require high-resolution space-based 
imaging and spectroscopic information for the individual ob- 
jects. 

4.4.3. Dark Matter halo properties 

The elongation in the total projected mass distribution seen in 
the weak lensing map and the high galaxy density associated 
with the group bullet suggest that the initial mass of the merg- 
ing group corresponds to a significant fraction of the main clus- 
ter. With the available ground-based WL data the group mass 
cannot be further constrained at the moment, but a reasonable 
assumption would be a bullet mass of the order of 10 14 Mq cor- 
responding approximately to a 1 :4 merger scenario. 

Hydrodynamical simulations to reproduce the ob- 
served properties of t he merging system IE 0657-56 by 
ISpringel & Farrail (120071) have shown that the time dependent 
spatial separation of the mass peak of the bullet group relative 
to the ICM and the shape of the contact discontinuity in the 



20 



R. Fassbender et al.: A pan-chromatic view of XMMU J 1230.3+ 1339 at z = 0.975 




Fig. 13. Zoom on the dynamically active North-Western core re- 
gion of the cluster. The 62"x55" LBT/LBC optical color com- 
posite displays the Northern part of core component (I) and a full 
view on region (II) identified in Fig. [10] The Chandra X-ray sur- 
face brightness contours (cyan) exhibit a pronounced cone-like 
structure with an approximate opening angle of 56 deg (yellow 
lines). Marginally extended 1.4 GHz radio emission is detected 
along the approximate symmetry axis of the configuration, dis- 
played by the green cross for the location, the white spatial er- 
ror ellipse, and the green ellipse indicating the observed source 
shape and extent. Character labels refer to candidate events of 
tidal stripping (t), galaxy-galaxy interaction (i), and a foreground 
galaxy (f). To enhance the contrast of low surface brightness fea- 
tures in the optical image, a smoothing with a 0.7" Gaussian 
kernel was applied to the data, and the black background was 
remapped to gray scale. 



X-ray surface brightness depend sensitively on the relative 
concentrations of the merging halos and their gas fractions. 
A first visual comparison of the merger configuration in 
XMMU J1230. 3+1339 to these simulations seem to be consis- 
tent with a low concentration parameter of c 2-3 for the main 
halo and a time snapshot within about lOOMyr of the closest 
core passage point. 

4.4.4. Observed extended radio emission 

An additional point of particular interest is the observation of 
extended radio emission originating from a location trailing the 
front of the infalling group very close to the approximate sym- 
metry axis of the merger configuration. In Fig.[T3] the location 
of the detected 1.4 GHz radio source (ID: J1230 1 6.1 + 133917) 
from the VLA FIRST survey (IWhite etalJll997t iBecker et all 
120031) is indicated by the green cross in the center of the image. 
This source with a signal-to-noise ratio of 7.2 is the only catalog 
entry within the field of the larger cluster environment shown in 
Fig. [10] The radio source is clearly extended with good signifi- 
cance along the North-South axis with a deconvolved FWHM of 
(5.7 ± 1.5)", but below the extent detection threshold along the 
orthogonal East- West axis. The reconstructed source shape is in- 
dicated by the green ellipse in Fig. [13] where the upper 1 cr limit 
of the extent along the minor axis of 1.5" was used. The white 



ellipse displays the positional error region of the source center 
at 90% confidence, which is just small enough to exclude any 
of the centers of the neighboring galaxies as origin for the radio 
emission at this confidence leveQ The 1.4 GHz radio source 
has an integrated flux of (1.69 + 0.19) mJy, which translates into 
a total radio power of P1440MHZ -8.3x10 24 WHz 1 at the cluster 
redshift under the assumpti on of a typical spectral in dex of -1 for 
the extended emission (e.g. lLaing & Peacocklll98dl) . 

4.4.5. What is the origin of the radio emission? 

The simplest explanation for the observed radio emission close 
to the core region of XMMUJ1230.3+1339 would be an ex- 
tended radio lobe associated with central AGN activity in one 
of the cores of the two red early-type galaxies partially en- 
com passed by the whi te positional error ellipse in Fig.Q~3] In 
fact. iDunn etail (1201 Oh have found extended radio emission in 
more than half of the early-type galaxies of a complete local 
sample of X-ray luminous E and SO galaxies. However, only 
the two most radio luminous objects in this sample (M87 & 
NGC 4696) would have comparable radio powers to the source 
in XMMUJ1230.3+1339, but these are associated with BCGs 
of local clusters. With the currently available radio data on the 
field, the association of the detected radio source with a nearby 
galaxy cannot be ruled out and is certainly a valid hypothesis. 

However, the fact that we are most likely witnessing an 
ongoing merger event close to cluster core passage suggests 
the consideration and discussion of an alternative scenario on 
the basis of turbulence-induced electron acceleration in the 
wake of infalling sub structure. This mech a nism of turbulence 
(re)acceleration (e g. iBrunetti et ail 1200 it iFujita et al.l 120031: 
iBrunetti et al.l2004T) is now believed to be the source of relativis- 
tic electrons associated with the giant radio halos observed in a 
growing number of m assive clusters with recent merger activ- 
ity su ch as Coma (e.g. iFeretti & Giovanninill2008t iFerrari et al.l 
120081 for reviews). The generation and time evolution of fluid 
turbulence behind a moving subcluster can be probed and visu- 
alized with numerical hydrodynamic simula tions (e.g. lTakizawal 
120051: iDolag et al.l [20081: IVazza et al.l [20091) . The overall geom- 
etry and involved length scales seen in such numerical studies 
are comparable to the observed configuration shown in Fig. [131 
implying that merger-generated turbulence at the location of the 
observed radio emission in XMMU J1230. 3+1339 is plausible 
and even expected. 

An intriguing scenario would be that we are actually witness- 
ing the injection point of re-accelerated ultra-relativistic elec- 
trons into the ICM as the source for the later development of 
a giant radio halo filling a good fraction of the cluster volume 
with low surface brightness radio emission. With the current ra- 
dio data, we can only perform some consistency checks to test 
the plausibility of t his alte rnative hypothesis. 

ICassano et al.l (120071) derived the expected scaling rela- 
tions of the re-acceleration model for the total emitted ra- 
dio power with total cluster mass and other properties, which 
is reflected in the observed strong empirical correlations. The 
latest radio p o wer c orrelation measurements are provided by 
Brunetti et al. (120091) for the Pi44omhz-L^' 1 relatior{3 yielding 
logtfWoMHz/WHz- 1 ) - 24.5 = (0.077 + 0.057) + (1.76 + 



15 The nearest previously known foreground galaxy (label f in Fig.l 1 3t 
SDSS J123015. 88+133920.4 has a spatial separation of 4.3" from the 
radio source, corresponding to a distance of about a factor of 4 outside 
the 90% positional error ellipse of the FIRST source. 

16 Using the BCES Bisector parameters. 
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0.16) [log(L^° /erg s ) — 45.4]. The combination of this correla- 
tion with the non-evolving Lt^ l 5Qij -Msm relation (Equ.© results 
in 
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Here f L f ,„„ =* 0.95 is a conversion factor from total lu- 

J tot— >500 

minosity to ^soo an< ^ /500->200 - 1-4 as defined in Sect.- 
13.1.11 Applying this relation to the measured radio power in 
XMMU J1230.3+ 1339, we obtain a mass estimate of M l 2 ^ GHz 
(9.5 + 3.8) x 1O 14 M (ME 13). The error is dominated by the 
systematic uncertainty in the redshift evolution of L^^-M 500 
(Sect. l3.l3l l. Other error contributions such as the uncertainties 
for the total radio flux (+12%), the K-correction from the as- 
sumed spectral index (+20%), and the used scaling relations are 
added in quadrature. 

This plausibility check for the hypothesis that the observed 
marginally extended radio emission is related to the origin 
of the radio halos in lower redshift clusters yields a mass in 
fairly good agreement with the best total mass estimate for 
XMMU J1230. 3+1339. Turning the argument around implies 
that for the inferred total cluster mass the observed total radio 
power is just a factor of a few higher than expected for the total 
radio halo emission, which might be expected close to the core 
passage of a subcluster. Hence, the scenario of the turbulence- 
induced electron re-acceleration model connected to the ongo- 
ing merging activity seems also to be a consistent alternative 
and should be further pursued. Considering that current detailed 
measur ements of diffuse rad io emission in clusters are limited to 
Z < 0.6 dBrunetti et alj200"9h . future low-frequency radio studies 
of the dynamically active core region of XMMU J1230. 3+1339 
have the potential to gain new insights into the origin and evolu- 
tion of radio halos. 



5. Discussion 

In the following section we discuss the dynamical state of 
the system and evaluate the different mass estimates and 
scaling relations used to derive them. We then compare 
XMMU J1230.3+1339 to other high-z galaxy clusters and focus 
on the connection to low redshift counterparts. 

5.1. Dynamical state of XMMU J1 230.3+1 339 

The system XMMU J1230. 3+1339 is a massive bona fide clus- 
ter with ongoing dynamical activity on all relevant scales. We 
will argue in the following that this cluster, observed at a look- 
back time of 7.6 Gyr, is likely a precursor of the most massive 
non-cool-core clusters in the local Universe, similar to the Coma 
system. 

In the optical, XMMU J1230.3+1339 (see Table© exhibits 
a galaxy population which is among the richest of all known 
z > 0.9 clusters. On the X-ray side, the ICM is found to be hot, 
with a high luminosity and a rather extended surface brightness 
profile, devoid of any detectable X-ray point source within Rsoo- 
Moreover, the high total mass of the system enables complemen- 
tary system diagnostics using weak lensing measurements. The 
cluster center of XMMU J1230.3+1339 is well defined, which 
indicates a virialized core region prior to the onset of the ob- 
served central merging event. The projected density maps reveal 
the galaxy and total light peaks, the center of the total mass, and 



the X-ray surface brightness peak all within a few arcseconds of 
each other (Fig.©. 



The location of the Brightest Cluster Galaxy at a projected 
distance of 18" or 140 kpc from this nominal center shows 
a significant spatial offset. Considering the dynamical friction 
timescale of a few Gigayears, spatial offsets of the BCG from 
the center of the main cluster halo are not surprising when ob- 
serving systems in the first half of cosmic time. However, in the 
local Universe the BCG location with respect to the X-ray cen- 
ter is kno wn to correlate strongly with the dynamical state of 
the cluster. ISanderson et al.l ((2009) have shown that a large pro- 
jected BCG offset (> 0.02 R500) implies on average a non-cool 
core host cluster with a shallow central ICM gas density profile, 
a larger than average gas mass fraction, and very low star for- 
mation and AGN activity associated with the BCG itself. These 
local (z ~ 0.2) correlations and expectations for clusters with 
significant BCG offsets are tentatively consistent with the prop- 
erties of XMMUJ1230.3+1339, with the exception of the gas 
mass fraction, which is found to be close to the average local 
value (see Sect. 13. 1 .5T >. 

The identified central merging event gives rise to a dynam- 
ically highly active central region of the cluster with various 
characteristic signatures as discussed in Sect. 14.41 Due to the 
quadratic dependence of thermal bremsstrahlung emissivity on 
the gas density, global ICM X-ray diagnostic measurements are 
expected to be influenced most by central group mergers. In the 
specific case of XMMU J1230.3+1339, the resulting flattening 
of the X-ray surface brightness within the wedge-like region in 
the wake of the infalling group is well visible in Fig. ©(top right 
panel). This feature is clearly a deviation from the underlying 
assumption of spherical symmetry and results effectively in a 
less reliable radial profile fit (Sect. l2.2^2l with an increased core 
radius r c and slope parameter /3 with respect to the pre-merger 
configuration. Additionally, significant deviations from an as- 
sumed hydrostatic equilibrium state (Sect. 13. It can be expected 
in the core region. The local ICM temperature could be lowered 
due to the cooler gas of the infalling group core, or temporarily 
boosted in the bow shock regions in conjunction with the X-ray 
luminosity. However, in the light of total cluster mass proxies 
of Sect. 15. 21 this merging event does not seem to affect the re- 
sults significantly. The interaction region is rather to be regarded 
as a localized perturbation with little influence on the still well 
defined cluster center and only mild effects on global quantities. 

Standard cool core diagnostics classify XMMU- 
J1230. 3+1339 in its current state as a Non-Cool-Core cluster 
(Sects. |2~2T21 & I3. 1 .41 ). albeit the relatively short central cooling 
time scale. Any possible previous onset for the development of a 
centrally peaked surface brightness profile was already flattened 
out by the group 'bullet' (component II). The foreseeable core 
impact of the BCG2 group (component III) will likely further 
disrupt the dense cluster center and hence erase the central 
conditions required for the progress of significant cooling. 
This scenario is consistent with the emerging picture from 
simulations on the different f ormation his t ories of NCC and 
CC clusters. The simulations of iBurns et al. I J20081) showed that 
Non-Cool-Core clusters are characterized by an increased mass 
accretion rate (^50% per Gyr) in early epochs (z ^ 0.8), during 
which nascent cool cores are destroyed by major merger events 
which also set the conditions to prevent cooling at later epochs. 
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Fig. 14. Comparison of different mass proxies for cluster 
XMMUJ1230.3+1339. The combined mass estimate Af*j* 
4.19*q™x10 14 Mq, shown on top with 1 cr uncertainties (dashed 
lines), is based on the combination of four quasi independent 
mass proxies (red points). The numbering along the vertical axis 
is equivalent to the mass estimate labels (ME 1-13) in the text. 
Statistical 1 cr uncertainties are shown by the black horizontal 
lines. The magnitude and direction of potential systematic off- 
sets are displayed by the slightly shifted blue bars extending the 
statistical error for illustration purposes. 



the cooling function for the soft X-ray band is practically in- 
sensitive to the ICM temperature (few per cent changes over 
the relevant Tx range). Similarly, the temperature based Z?200 
was chosen as cluster boundary for the weak lensing measure- 
ments, which was regarded as more robust compared to the self- 
consistent approach based on the current model constraints from 
the ground-based weak lensing signal with larger uncertainties. 
We choose the NFW model estimate (ME 12) as WL input mass 
proxy for the MSj? determination. As fourth reference mass, the 
'optimal optical' mass tracer is considered (ME 9). 

Besides the statistical uncertainties arising from measure- 
ment errors and the scatter in the used scaling relations, we also 
attempted to approximately quantify the influence of potential 
systematic offsets on the mass estimators. These latter uncertain- 
ties were due to e.g. smaller analysis apertures compared to the 
used scaling relations (Tx), the resulting influence on the derived 
cluster radii (WL), significant uncertainties in the redshift evo- 
lution (Lx, Lbcg, optical apertures), and technical ambiguities 
arising from the adaptation of local measurement recipes to high 
redshift (Nwo, £200) ■ These technical uncertainties are in prin- 
ciple independent from the statistical errors, however, their (un- 
known) statistical occurrence may not follow a Gaussian prob- 
ability distribution and are hence difficult to properly account 
for. As a first order approximation, we assume Gaussian prop- 
erties of the systematic technical uncertainties and add them in 
quadrature to the statistical error intervals to obtain mass prob- 
ability distributions for the individual proxies. The foufl quasi- 
independent mass reference proxies are then combined with a 
maximum-likelihood approach, resulting in a final cluster mass 
estimate of = 4.19+° ™ x 1O 14 M with specified 1 cr un- 
certainties. This currently most likely total mass interval for 
XMMU J1230. 3+1339 is shown in the mass proxy comparison 
plot in Fig.[14](vertical dashed lines), together with the input ref- 
erence proxies displayed by the red filled circles, and all other 
mass estimates (black circles) with their corresponding uncer- 
tainties. 



5.2. Comparison of mass estimates and scaling relations 

In Sects. [3] & I4.4I we have derived a total of 13 different mass 
proxies for XMMU J1230.3+ 1339, which will be compared and 
discussed in the following. 

5.2.1. Combined mass estimate 

We have applied different and independent observational tech- 
niques to study and characterize the cluster based on (i) opti- 
cal imaging data, (ii) X-ray observations, (iii) weak gravitational 
lensing, and (iv) spectroscopic information for a first glimpse on 
the system dynamics. The objective in this section is to combine 
several suitable measurements in order to obtain a (preliminary) 
best total mass estimate as a benchmark value and starting 
point for a comparison of the different high-z mass proxies. 

We start by selecting four reliable, quasi-independent mass 
estimates from different observational techniques as input for 
a combined mass estimate. On the X-ray side, the concordant 
Chandra and XMM-Newton ICM temperature estimate provides 
tight mass constraints based on the Tx-M (ME 1) scaling rela- 
tion, which also encompasses a reliable redshift evolution pre- 
diction. The Chandra determined /3-model SB structure parame- 
ters were the basis for the gas mass measurements of the system 
and the corresponding M|™ estimate (ME 4). This constitutes a 
second reliable and quasi-independent mass measurement, since 



5.2.2. Multi-wavelength mass proxies and challenges at 
high-z 

The combined cluster mass proxy - 4.2 X 10 14 Mq with 
l cr uncertainties of about 19% is currently the best reference 
value for XMMU J1230. 3+1339 and is among the most accurate 
cluster masses measured to date at z > 0.9. 

The X-ray derived mass estimates (ME 1-5) provide a con- 
sistent picture of the cluster with mostly well constrained mass 
intervals based on the medium deep Chandra and XMM-Newton 
data. The systematic uncertainties on Tx (ME 1) arise from the 
unconstrained temperature structure of the ICM and the possi- 
ble influence of the central merging event. The constraints from 
the hydrostatic estimate (ME 2) are mostly limited by the un- 
certainties in the /^-parameter, which is certainly influenced by 
the fly-through group bullet in addition to possible deviations 
from the hydrostatic assumption in the central regions. The un- 
certainties of the gas mass (ME 4) and Yx (ME 5) estimates are 
mostly limited by the depth of the available X-ray data. The clus- 
ter property with the smallest relative uncertainty (see TableO 
is the X-ray luminosity Lx, which translates into small statisti- 
cal mass estimate errors based on the latest (local) Lx-M scal- 
ing relation. However, a significant uncertainty is added because 



The velocity dispersion-based dynamical mass estimate (ME 10) 
with the current uncertainty range has practically no statistical weight, 
its inclusion for the combined mass would only have a per cent effect. 
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of the limited knowledge on the high-z evolution of the related 
L-T scaling relation (see Sect. l3.1.3l l. For the very good concor- 
dance of the Lx based mass estimate with M™, the Lx-T rela- 
tion was assumed to be non-evolving, which is different from the 
self-similar evolution expectations. However, such a breaking of 
self-similarity for the Lx-M scaling in terms of a slower redshift 
evolution and a steeper slope of the relation is in qualitatively 
good agreement with the latest results from the Millennium Gas 
Simulations a fter the inclusion of non-gravitational (pre-)heating 
mechanisms (IStanek et al.l l2010h . Once an improved and self- 
consistent calibration of the Lx scaling relations at z > 0.9 is 
available, the X-ray luminosity appears to have the potential to 
be a good high-z mass proxy at relatively low observational cost. 

Weak lensing mass estimates are potentially a powerful tool 
to provide independent low -bias measurements to make progress 
in the calibration of mass-observable relations. However, the 
weak lensing analysis of XMMU J1230.3+1339 in Paper II is 
at the current feasibility limit for ground-based imaging obser- 
vations, implying fairly large statistical uncertainties due to the 
limited lensing signal strength. M™' is on the lower mass end of 
the uncertainty range but still consistent within the errors with 
the WL mass proxies (ME 1 1 & 12). Significant progress and a 
self-consistent parameter treatment in the WL analysis will re- 
quire the availability of deep spaced-based imaging data. 

The dynamical cluster mass estimate based on the measured 
radial velocity dispersion o> (ME 10) provides the weakest sta- 
tistical constraint for XMMU J1230.3+1339, due to the small 
dynamic range of <x,. and the relatively large uncertainties related 
to the limited number of available redshifts. The situation of hav- 
ing about a dozen spectroscopic cluster members for a first or- 
der evaluation of the velocity dispersion is typical under survey 
conditions with single mask spectroscopic follow-up of newly 
discovered high-z systems. More stringent mass constraints re- 
quire large, dedicated follow-up programs to ideally increase the 
number of cluster redshifts by an order of magnitude. Although 
hardly feasible in the general case for z > 0.9 clusters due to the 
magnitude limit of the available spectrographs, the high richness 
of XMMU J1230.3+1339 might realistically enable the identi- 
fication of about 100 spectroscopic cluster members. However, 
biases arising from infalling structures, which could lead to pos- 
itive and negative offsets depending on the line-of-sight angle, 
are likely to limit the achievable accuracy of dynamical mass 
estimates at high redshifts. 

Mass proxies based on optical/NIR imaging data have the 
advantage of being observationally cheap and usually readily 
available for distant cluster searches or studies. Although sig- 
nificant progress has been made in the calibration of optical 
scaling relations at lower redshifts, unbiased mass proxies at 
z > 0.9 to better than 40-50% seem to require significant fur- 
ther work to overcome technical and evolutionary challenges. 
We have obtained a consistent picture that the optical galaxy 
population of XMMU J1230.3+1339 is rather remarkable for 
a high-z cluster, both in terms of total galaxy numbers (A^oo, 
R nch ) and total luminosity (L200, L^L). Interestingly, the optical 
BCG luminosity (ME 7) yields a mass estimate which is fully 
consistent with M™, followed by the suggested higher system 
masses of the A^oo proxy (ME 6), the proposed optimal opti- 
cal mass tracer (ME 9), and the L200 luminosity (ME 8). In the 
ca se of XMMU J123 0. 3+1339, the optimal optical mass tracer 
of iReves et all d2008l) does not seem to improve the individual 
proxies based on A^oo and Lbcg- 

The use of luminosity -based mass proxies (e.g. L200) at high- 
z requires as input rather large evolution and K-correction fac- 
tors to relate them to the scaling relations calibrated at lower 



redshifts. Counting galaxies along the typically well-defined red- 
sequence of early-type galaxies (e.g. A^oo) seems to be a more 
straightforward approach. However, even this apparently sim- 
ple task gets increasingly complicated at high-z, where various 
subtleties and potential biases are to be taken into account if a 
somewhat reliable mass estimate is to be achieved. Even under 
the assumption that a distant cluster is sufficiently rich to exhibit 
an easily identifiable red-sequence, three main critical techni- 
cal aspects remain and may impose significant mass biases: (i) 
finding the correct RS slope, (ii) the definition of appropriate 
color cuts, and (iii) the determination of a self-consistent mea- 
surement aperture. All three items are in principle functions of 
cluster redshift, and the first two depend moreover on the ap- 
plied filter combination and the depth and quality of the data. As 
examples of redshift dependent biases, we can consider galaxy 
evolution effects in clusters, which lead to an observed increase 
with lookback time of the fraction of objects bluer than the red- 
sequence and also an increasing relative effective age (and hence 
color) difference of the average stellar population s as a func- 
tion of cluster-centric radius (e.g. lRosati et alj|2009t) . For a fixed 
color cut definition, a redshift dependent positive or negative 
richness and mass bias could thus be imposed through an in- 
creasing contamination of bluer non-passive galaxies in the red 
ridgeline selection, or by excluding a growing fraction of slightly 
bluer passive galaxies on the cluster outskirts. Similar redshift 
dependent mass biases may arise through a non-consistent scal- 
ing of the measurement apertures. 

An alternative approach applicable in the absence of a robust 
estimate for the cluster radius 7?2oo or for direct cluster-to-cluster 
comparisons is to measure relevant cluster quantities inside a 
fixed physical projected radius of e.g. 1 hzl Mpc. The hereby im- 
posed positive mass bias with redshift might in many cases be an 
acceptable trade-off for otherwise large possible uncertainties in 
the cluster radius. In the case of XMMU J1230.3+1339, the ICM 
temperature determined /?200 is within 2% of 1 Mpc, imply- 
ing effectively only percent level differences for most of the total 
mass estimates (ME 1-5, 7, 10-13), when re-scaled to the fixed 
physical aperture. However, for the optical mass proxies based 
on richness and total luminosity (ME 6, 8, 9) a significant re- 
duction of the total mass estimate of the order of 30% applies 
for the fixed radius, which would move these proxies in much 
better agreement with the current best estimate M™. 

5.3. Comparison to similar X-ray clusters atz~l 

The growing number of known z ~ 1 galaxy clusters with avail- 
able multi-wavelength data will soon allow comprehensive pop- 
ulation studies from the group regime to the most massive sys- 
tems at this important cosmic epoch. A point of particular inter- 
est will be the intrinsic cluster-to-cluster variations at a given 
epoch for systems with comparable masses, e.g. with respect 
to the thermodynamic properties of the ICM, the formation his- 
tory of the galaxy population, or the structure of the Dark Matter 
halo. This intrinsic physical scatter sets on one hand the ultimate 
limit of accuracy for the calibration of high-z scaling relations 
and reflects on the other hand the different conditions during 
cluster formation, e.g. the collapse epoch, the influence of the 
LSS environment, or AGN feedback. 

Within the XDCP survey, the cluster XMMU J 1229.5 +0151 
at z = 0.975 (ISantos et al.ll2009l) is almost the 'twin brother' 
of XMMUJ1230.3+1339. Besides the concordant redshift, 
XMMU J1229.5+0151 features a comparable X-ray luminosity, 
a similar ICM temperature within the errors (Tx — 6.4 keV), 
and a consistent velocity dispersion (cr r = 683 km/s). The clus- 
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ter is also optically rich with a well populated red sequence and 
two bright galaxies in the center. A third bright galaxy of com- 
parable magnitude and mass is located on the outskirts of the 
cluster, at a projected distance similar to the BCG2 group in 
XMMUJ1230.3+1339. Overall, the basic global properties of 
the two clusters seem to be well consistent within the observa- 
tional uncertainties. 

The two most distant clusters from the Wide A ngle 
ROSAT Pointed Survey (WARPS) dPerlman et all 120021) are 
further appropriate candidates for a direct comparison with 
XMMU J1230. 3+1339 as they are similar in mass and redshift. 
The most distant WARPS cluster, CI 1415.1+3612 at z = 1.03, 
has an ICM temperature of Tx — 5.7 keV, a luminosity of 
L^°2 00 - 10.4 x 10 44 erg/s, a hydrostatic total mass estimate of 

M 2 oo - 3.8 x 10 14 M© dMaughan et al.l 120061). a nd a velocity 
dispersion of <x r - 810km/s dHuang et al.l |2009|) . The second 
system CI 1429.0+4241 at z = 0.92 is characterized by the pa- 
rameters T x 6.2 keV, L^ 00 ^ 9.6 x 10 44 erg/s, and M 20 o - 

4.5 x 10 14 M Q . With core radii of r c ~ 100 kpc (for /? = 2/3) both 
systems are more compact in their ICM structure compared to 
XMMU J1230.3+1339 and their bolometric X-ray luminosities 
are 40-50% higher, while the total mass and temperature esti- 
mates are similar within the uncertainties. 



5.4. Appearance atz-0 

We now address the question of the expected system properties 
of XMMUJ1230.3+1339 after an additional 7.6Gyr of cluster 
evolution (i.e. z = 0) and compare these to the characteristics of 
the local Coma cluster (z- 0.023). 

A lower limit for the total system mass at z = 0, without 
considering further accretion, can be obtained from the self- 
similar evolution model. Due to the decreasing reference den- 
sity p a (z) with time, the system radius is expected to evolve 
as /?20o(z = 0) - E 2/i (z) ■ /?20o(z), which translates approxi- 
mately linearly into the total cluster mass Moodiz -0)- E 2 ^{z) • 
Moooiz) - 1.4 • M 2 oo(0.975) (see e.g. Equ.©. However, once in- 
cluding the continuing mass accretion the cluster growth rate 
since z - 1 is much faster. The simulations of iBurns et al. I J2008T) 
predict average mass growth rates of ^5-10% per Gyr at late 
times and ^30-50% per Gyr during active mass assembly pe- 
riods around z ~ 1 . As a benchmark, the time evolution of the 
largest local ~10 15 Mp hal o in the Millennium-II simulation 
dBovlan-Kolchin et al.ll2009l) exhibits a total mass increase since 
z ~ 1 by more than a factor of 5. A realistic prediction for the 
evolved total cluster mass of XMMU J1230.3+1339 at the cur- 
rent epoch is hence M 2 oo(z = 0) ~ (1 .5-2) x 10 15 Mq, i.e. ranked 
among the most massive local clusters and likely even higher in 
mass compared to the Coma system with 10 15 Mq. 

From the X-ray perspective, Coma is a prototype Non-Cool- 
Core cluster featuring a shallow surf ace brightness pr ofile with a 
large core radius of r c ~ 300 kpc (e.g. lBriel et alj 1992h . A similar 
X-ray appearance can be expected for XMMU J1230. 3+1339, 
which shows already in the current configuration a very extended 
NCC SB distribution. The core impact of the BCG2 group and 
subsequently accreted matter (Sect. 15. it will likely further dis- 
perse the remaining central ICM density peak within the next 
Gigayear (Figs. [8] [9]) resulting in an even more flattened X-ray 
SB profile for the expected new equilibrium configuration a few 
Gyr later. 

The upcoming merging with the BCG2 group will lead to 
another close resemblance to Coma in the optical, namely the 
existence of two central, almost equally massive brightest galax- 



Table 3. Multi-wavelength 
XMMUJ 1230+1339. 



properties of the cluster 



Property Value 



Unit 



M%,[10"Mq] 



RA 
DEC 

N H 

z 

lyjbcst 



Rx,200 
Rx,500 
Rx.2500 



12:30:16.9 
13:39:04 
2.62 x 10 20 
0.975 ± 0.002 
4.19^ x 10 14 

Ton** 

670^ 
297 + ?° 



Up 







kpc 
kpc 
kpc 



Tx,2500 


c an+0.71 
-0.62 


keV 


3 - 61 -0.68 


l x,n" 


f. n+i 6 


Ke v 




■^X.2500 


43+ - 24 






f0.5-2keV 
X.500 


(5.14 + 0.54) x 10~ 14 


erg s cm -2 




t 0.5-2 kcV 
X.500 


(1.92 ± 0.20) x 10 44 


ers s' 




t 0.5-2 kcV 
^X,[0.15-l]500 


(1.53 + 0.21) x 10 44 


erg s - ' 




T bol 

L X,500 


(6.50 ± 0.68) x 10 44 


erg s - ' 


4 1 ±0.6 (+1.0) 


Z2500 


43+0.24 


z 




fix 


0.84 ± 0.47 




4.9 ±2 - 8 


IV 


215 ± 110 


kpc 




C SB 


0.07 






"e.O 


(2.3 + 1.4) x 10~ 2 


cm -3 




tcool,20kpc 


4.0 ± 2.4 


Gyr 




''sound 


1160 


km/s 




M gaS- 500 


(3.0 ± 0.9) x 10 13 


M Q 


42+1.5 


Yx,500 


L6± o.5 (+ o.2) x 10 i4 


M keV 


36 ±l.l(+0.4) 


fgas,500 


0.100 + 0.035 




R nch 


2 






N200 


57 




g 3 ±0.7 (±2.0) 


^200 


(9.2 + 2.3) x 10 12 


Lq 


lj j±2.8(±3.7) 


j tot.z 
.500 


(1.06 ± 0.16) x 10 13 


L© 




R-z 


2.05+0.10 


mag 





m BCG 

JV1 BCG 
^center 
U BCG 
CT,. 

M/L Z ' 
Optopt 



19.90 + 0.02 
-24.95 ± 0.03 
140 

658 ± 277 
46.7 ±11.3 



mag 
kpc 
km/s 
Mo/L© 



< 4 +0.7 (+1.0), 
^' -0.7 (-2.2) > 



2.8I| 2 



8.6 



WLoj ls 
WL NFW 



1271 ±275 



km/s 



1440MHz 



8.3 x 10 2 ' 



WHz~ 



+1.3 (+2.8) 

1.3 (-3.9) 
7 fi +3.3(+3.0) 

'•"-3.3 
f- q+3.3(+2.2) 
"• 7 -3 3 
-95^ 



ies. The core impact velocity of the BCG2 group with the main 
halo is too high (see Sect. l4.41 > to allow efficient merging with 
other central galaxies, i.e. BCG2 is expected to retain its en- 
tity as a massive galaxy. After relaxation within the main clus- 
ter halo and several Gyr of dynamical friction at work, BCG2 
can be expected to settle close to the cluster center at low red- 
shifts and form a pair of two dominant cluster galaxies together 
with the nominal BCG In terms of the general optical richness, 
XMMU J1230. 3+1339 is already comparable to Coma (richness 
class R=2) when considering the accessible bright end of the 
galaxy population. With the continuing accretion of galaxies 
from the surrounding filaments, the new system is likely to sur- 
pass Coma in optical richness by z-Q. 

Coma is also the prototype for clusters featuring a diffuse 
radio halo. In Sect. 14.41 we discussed the possibility that such a 
radio halo originates from merging events such as the one ob- 
served in XMMU J1230.3+1339. If this scenario prevails, the 
evolved cluster might additionally show similarity to Coma in 
the radio regime at certain epochs. 
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6. Summary and conclusions 

We have presented a multi-wavelength analysis and characteri- 
zation of the galaxy cluster XMMU J1230.3+1339 at z = 0.975, 
which is summarized in the following. An overview of the sys- 
tem parameters and the derived mass estimates is provided in 
Table[3] 

1 . The galaxy cluster was discovered within the XMM-Newton 
Distant Cluster Project (XDCP) as a bright, extended X-ray 
source in an archival XMM-Newton observation targeting 
the Virgo galaxy NGC 4477. 

2. We have characterized the X-ray properties of the system 
with a joint analysis of the XMM-Newton data and addi- 
tional archival Chandra observations finding a bolometric 
luminosity of L]^ m = 6.50 x lO^ergs" 1 (± 10%), an ICM 
temperature of 7x,2500 — 5.30 keV (+ 13%), a gas mass of 
M gas> 5oo - 3.0 x 10 13 M Q (+30%), a gas mass fraction of 
fgas,500-0.10 (+ 35%) consistent with local values, and a Yx 
parameter of Fx,500 - 1 -6 x 10 14 M keV (+ 38%). 

3. The Chandra analysis of the X-ray surface brightness dis- 
tribution yielded an approximate core radius of 215 kpc, a 
f3 profile parameter of 0.84, and an overall classification as 
a Non-Cool-Core cluster. The very central part of the SB 
profile is better represented by a double /3-model with an 
additional more compact radial component, from which we 
derive a central electron density of n e ,o - 2.3 x 10~ 2 cm~ 3 
(+ 60%) and a corresponding cooling time of about 4Gyr. 

4. The cluster features a very rich galaxy population (Abell 
richness R=2, A^oo - 57) with a prominent red-sequence in 
the color-magnitude-diagram and a total z'-band luminosity 

of^ p -l.06xl0 13 L o (+l5%). 

5. The Brightest Cluster Galaxy is located at a projected dis- 
tance of 140 kpc from the nominal cluster center and has a 
luminosity of L Z BCC 6.1x10" Lq, two magnitudes brighter 
than the characteristic L* luminosity at this redshift. 

6. Based on 13 secure spectroscopic cluster members, we find a 
radial velocity dispersion estimate of cr r 658 km/s (+42%). 

7. As presented i n deta il in the accompanying Paper II of 
iLerchster et alJ d2010h . we detect a weak gravitational lens- 
ing signal for the cluster on the 3.5 cr significance level, 
which currently sets the redshift limit for the extraction of 
a weak lensing signal based on ground imaging. 

8. We have performed a multi-wavelength analysis of the clus- 
ter structure based on projected density maps of (i) red galax- 
ies, (ii) the X-ray surface brightness, (iii) the total z'-band 
light, and (iv) the weak lensing mass distribution. We find a 
well defined concordant cluster center and a main elongation 
axis along the SE-NW direction present in all maps. 

9. Several cluster-associated sub-components were identified 
on different cluster-centric distance scales ranging from a 
merging event in the central region to two infalling groups 
on the cluster outskirts and two galaxy filaments extending 
beyond the nominal cluster radius. 

10. One of the groups has recently entered the R^q radius of 
the cluster along a radial infall orbit. The central domi- 
nant galaxy of this group has also BCG properties and is 
only marginally fainter than the nominal Brightest Cluster 
Galaxy. 

11. We find evidence for a central fly-through group 'bullet' 
close to cluster core passage from a wedge-like X-ray sur- 
face brightness feature, an associated off-center galaxy den- 
sity and luminosity peak with the front end aligned to the 
X-ray surface brightness jump, indications for tidally dis- 



rupted galaxies, blue-shifted redshifts in the cluster-centric 
reference frame, and a matching elongation in the weak lens- 
ing mass map. 

12. Marginally extended 1.4 GHz radio emission located in the 
wake region of the central merging group along the axis 
of the apparent 'bullet' trajectory is detected by the VLA 
FIRST survey. The origin of this radio emission could be as- 
sociated to an extended radio jet connected to the center of 
one of the nearby cluster galaxies, or possibly the direct ob- 
servable result of ultra-relativistic particle acceleration in the 
turbulent wake region of the merger, as predicted by simula- 
tions. A cross-check with the scaling relation for the total 
power in lower-z radio halos with the observed radio source 
in the cluster is reasonably consistent with the possibility that 
they share a related physical origin. 

13. We have attempted a tentative reconstruction of the 3D- 
accretion geometry of the cluster environment that connects 
the identified associated components with available redshift 
information to the main halo. 

14. We have derived 13 proxies for the total cluster mass from 
the different wavelength regimes. The combination of four 
reliable quasi-independent mass estimates yielded a joint 
best estimate for the total cluster mass of Mfe? 4.19 x 
10 14 Mq (± 19%). All X-ray derived mass estimates are in 
good agreement with this preferred mass interval, the weak 
lensing results are consistent within the currently larger sta- 
tistical uncertainties, and the tested optical mass proxies tend 
to suggest significantly higher total system masses, albeit 
with large systematic uncertainties. 

15. We argued that the cluster will likely evolve into a Coma- 
like system with expected similarities in the ICM structure, 
the optical appearance with two dominant central galaxies, 
and possibly even surpass Coma in terms of galaxy richness 
and total mass. 

The newly identified galaxy cluster XMMU J1230.3+1339 is an 
excellent test laboratory at z ~ 1 for detailed observational stud- 
ies to address various fundamental questions related to the as- 
sembly phase of massive systems and their Dark Matter, ICM, 
and galaxy components. Recent follow-up campaigns of the 
cluster environment include additional optical spectroscopy at 
VLT/FORS 2, deeper X-ray coverage with Chandra, and mm 
observations with APEX-SZ. Upcoming studies will e.g. aim to 
derive a more detailed picture on the interaction processes of 
the different components related to the central merging event or 
to characterize the galaxy populations as a function of cluster- 
centric distance in different sub-components and the large-scale 
structure filaments. 
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